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ABSTRACT 


This  report  covers  the  results  of  limiting 
guidance  accuracy  studies  for  nine  proposed  aided 
inertial  guidance  techniques.  In  each  case  an  opti¬ 
mum  information  handling  scheme  was  assumed.  This 
provides  a  consistent  framework  for  comparing  alter¬ 
native  schemes  and  determining  the  desired  lower 
limit  on  errors  for  each  technique.  The  report  also 
develops  the  sensitivities  of  missile  guidance  errors 
to  sensor  errors  and  initial  alignment  errors  for  three 
unaided  guidance  techniques  —  inertial  guidance  and 
heading  command  guidance  with  and  without  wind  com¬ 
pensation.  The  sensitivity  tables  provided  allow  rapid 
comparisons  of  different  inertial  sensors  and  align¬ 
ment  techniques  which  may  be  suggested  for  use  ih 
proposed  guidance  systems. 
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1.  INTRODUCTION 


This  report  concerns  a  continuing  study  of  various  air-to-sur- 
face  missile  guidance  techniques.  It  presents  the  results  of  limiting  ac¬ 
curacy  studies  for  various  inertial  guidance  techniques.  It  also  presents 
the  sensitivity  analysis  of  missile  guidance  errors  to  sensor  errors  and 
initial  alignment  errors  for  heading  command  guidance,  position  reference 
(inertial)  and  wind  compensated  heading  command  guidance. 


1.1  OVERVIEW 

Accuracy  Limits  for  Aided  Inertial  Guidance  Techniques  — 
Every  new  guidance  technique,  whether  based  on  a  new  advance  in  sensor 
design  or  a  new  data  blending  scheme,  has  associated  with  it  basic  per¬ 
formance  limitations.  It  is  important  that  these  limitations  be  evaluated 
in  each  case  and  that  decisions  not  be  based  on  oversimplified  analyses 
which,  for  example,  treat  only  the  reductions  in  error  contributions 
which  previously  limited  performance.  Calculation  of  lower  bounds  on 
the  attainable  accuracy  of  specific  guidance  techniques  requires  identi¬ 
fication  of  fundamental  error  sources  and  determination  of  their  individual 
contributions  to  system  error.  These  calculations  may  be  simplified  by 
assuming  that  the  system  under  study  employs  an  optimum  information 
handling  scheme.  This  assumption  provides  a  consistent  framework 
against  which  alternative  suboptimal  schemes  can  be  compared.  It  also 
provides,  for  each  competing  guidance  technique,  the  desired  lower  limit 
values.  This  is  the  viewpoint  adapted  '  the  limiting  accuracy  studies 
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reported  in  this  document.  The  limiting  accuracies  of  nine  proposed 
inertial  guidance  techniques  are  discussed.  These  systems  cover  a  wide 
range  of  possible  system  designs.  The  limiting  factors  considered  can 
be  grouped  into  the  following  areas:  geometric  dilution,  propagation 
uncertainties,  geophysical  uncertainties,  electronic  clock  accuracy, 
and  other  external  aid  fix  errors.  The  effect  of  these  factors  on  the 
guidance  techniques  were  considered  and  detailed  numerical  results  are 
presented. 

Inertial  Guidance  Sensitivity  Study  —  Sensitivity  studies  were 
performed  relating  air -to -surface  missile  guidance  errors  to  sensor 
errors  and  initial  alignment  errors  for  five  representative  trajectories. 
The  guidance  techniques  considered  are: 

•  Heading  Command 

•  Position  Reference  (Pure  Inertial) 

•  Wind  Compensated  Heading  Command 

Results  are  presented  in  the  form  of  sensitivity  tables  which  allow  the 
rapid  evaluation  of  guidance  errors  for  particular  trajectories,  based 
on  inertial  sensor  errors  and  alignment  errors.  These  tables  are  valu¬ 
able  in  performing  accuracy  tradeoff  studies  of  inertial  sensors  and 
alignment  techniques  for  use  in  proposed  guidance  systems. 

1.  2  FORMAT  OF  THE  REPORT 

Chapter  2  documents  the  limiting  accuracy  studies  of  the  aided 
inertial  guidance  techniques.  The  limiting  factors  are  discussed  and  the 
results  of  the  analyses  of  the  nine  techniques  are  presented.  This  chapter 
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is  augmented  by  Volume  II  (classified  SECRET)  of  this  report,  which 
provides  details  on  satellite/inertial  and  radar  correlation/inertial 
techniques.  Chapter  3  describes  the  guidance  accuracy  sensitivity  anal¬ 
yses  performed.  The  three  guidance  techniques  and  five  trajectories 
considered  are  described.  The  sensitivity  tables  for  each  of  the  tra¬ 
jectory/guidance  technique  combinations  are  given  with  sample  calcu¬ 
lations  demonstrating  their  use.  Chapter  4  summarizes  the  main  findings 
of  the  work  described  in  this  document  and  relates  the  value  of  the  anal¬ 
yses  to  current  and  future  air-to-surface  missile  guidance  problems. 
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2.  ACCURACY  LIMITS  FOR  AIDED  INERTIAL 

GUIDANCE  TECHNIQUES 

» 

This  chapter  concerns  the  calculation  of  the  limiting  accuracy 
for  nine  proposed  inertial  guidance  techniques  for  air -to -surface  missiles. 
These  accuracy  limits  provide  both  a  lower  bound  on  the  attainable  ac¬ 
curacy  of  the  specific  guidance  techniques  and  a  consistent  baseline  for 
comparison  of  alternative  schemes.  The  factors  considered  to  limit  the 
guidance  accuracy  were  geometric  dilution  propagation  uncertainties, 
geophysical  uncertainties,  electronic  clock  accuracy,  and  other  external 
aid  fix  errors.  These  factors  are  discussed  in  detail  in  the  next  section. 
The  influence  of  these  factors  on  the  nine  inertial  techniques  listed  below 
is  also  contained  in  this  chapter. 

•  Pure  Inertial  Guidance 

•  LORAN/lnertial  Guidance 

Direct  Ranging  LORAN/'bertial  Guidance 

•  OMEGA/lnertial  Guidance 

•  DME/lnertial  Guidance 

•  Satellite/lnertial  Guidance* 

•  Radar  Correlation/lnert'.al  Guidance* 

•  Optical  Correlation/lnertial  Guidance 

•  Doppler/lnertial  Guidance 

$ 

Details  are  presented  in  the  Volume  II  of  this  report  (Ref.  7). 
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/  R  *  Tru*  Receiver  location 


'  1  *  Indicated  Receiver  location 

£f,0f*  *n  Determining  IOP  1  ond 
LOP  2  Respectively 
d  »  Radiol  Receiver  location  Error 
0 »  Angle  o(  Crossing  Between  the 
Two  LOPS 

Figure  2. 1-1  Effect  of  LOP  Errors  on  Receiver  Location  Error 


Assuming  Cj  and  are  zero-mean  random  variables  whose  statistics 
are  described  by  the  following  quantities: 

*  standard  deviation  of  the  error,  in  the 
measurement  of  LOP  1 

cr„  =  standard  deviation  of  the  error,  c_,  in  the 
measurement  of  LOP  2 

£  > 

p  »  the  correlation  coefficient  between  the  mea¬ 
surement  errors  of  LO.y  1  and  LOP  2 

and  taking  the  expected  value  of  the  terms  in  Eq.  (2. 1-3),  the  square  of 

2 

the  standard  deviation  of  the  radial  error  in  receiver  location,  o^,  is 
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= 


1/22  \ 
~2~  lal  +a2  +  2Paicr2cos  6J 


sin  0 


(2,1-4) 


"  sin?  (“l +ff2  +2(wi°2cos  e) 


(2. 1-5) 


As  can  be  seen  from  Eq.  (2. 1-5)  the  standard  deviation  of  the 
receiver  location  error  is  highly  sensitive  to  the  angle  between  the  LOP’s. 
This  geometry  dependence,  which  has  the  effect  of  always  raising  cr^, 
when  0-ff/2  is  not  a  multiple  Of  rr,  is  known  as  geometric  dilution. 

Two  major  types  of  radio  navigation  systems  are  considered 
here.  One  generates  hyperbolic  lines  of  position  as  in  Fig.  2. 1-2  by 
measuring  phase  or  time  differences;  examples  are  OMEGA  or  conven¬ 
tional  LORAN  systems,  tt  can  be  shown  that  the  angle  of  crossing  of  the 
LOP’s  is 


4>,+«.2 


(2. 1-6) 


for  the  three -station  configuration  shown  in  Fig.  2. 1-2.  The  second  type 
of  radio  navigation  system  considered  generates  circular  lines  of  position, 
as  in  Fig.  2. 1-3,  by  providing  direct  measurement  of  range  to  the  receiver 
(one  possibility  is  that  the  measurement  is  in  the  form  of  time  for  the  signal 
to  arrive  from  the  station,  which  can  be  converted  to  range  simply  by  multi¬ 
plying  the  propagation  time  by  the  velocity  of  light)  from  one  station,  as 
with  Direct  Ranging  LORAN  (see  Ref.  1).  For  this  type  of  system  the  angle 
of  crossing  of  the  LOP’s  is 


9  *  (ff  -  i p) 


(2. 1-7) 


where  <p  is  defined  in  Fig.  (2. 1-3). 
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By  combining  Eqs.  (2. 1-5)  and  (2. 1-6)  it  can  be  shown  that 
for  a  hyperbolic  system  the  standard  deviation  of  the  radial  error  in 
receiver  location,  o^,  is  given  by 


+02+2^0^2  cos 


1/2 


(2. 1-8) 


where  <p^  and  are  determined  by  the  fix  geometry  (Fig.  2. 1-2).  Com¬ 
bining  Eqs.  (2. 1-5)  and  (2. 1-7)  gives  the  rms  radial  error  for  a  circular 
system  as 


sin<p 


(°i 


+  0^2  cos 


*) 


1/2 


(2. 1-9) 


where  <p  is  defined  by  the  fix  geometry  (Fig.  2. 1-3). 

2. 1. 2  Propagation  Uncertainties . 

Radio  waves  can  be  considered  to  propagate  over  the  earth 
in  either  a  ground  wave  and/ or  skywave  mode.  Ground  waves  propagate 
by  following  the  curvature  of  the  earth  and  are  generally  high  frequency 
waves.  This  mode  of  propagation  is  characteristic  of  such  systems  as 
LORAN.  Sky  waves  propagate  by  reflecting  in  a  wave  guide  formed  by 
the  ionosphere  and  the  earth’s  surface.  These  waves  are  generally 
low  frequence,  characteristic  of  such  systems  as  OMEGA.  The  two 
propagation  modes  are  illustrated  in  Fig.  2. 1-4.  Since  the  ground 
waves  are  propagated  directly  to  the  receiver  the  time  between 
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transmission  and  reception,  T,  can  be  related  to  the  range  (ground 
distance  or  slant  range),  R,  between  the  transmitter  and  receiver  by 
Eq.  (2.1-10) 

T  =  ?.  (2.1-10) 

c 

where  c  is  the  velocity  of  light.  The  latter  is  equal  to  the  nominal 
velocity  of  light  in  the  atmosphere,  cQ,  and  the  variation  in  the  velocity, 
6c,  due  to  uncertainties  in  atmospheric  conductivity  caused  by  atmo¬ 
spheric  disturbances,  water  vapor  concentration,  etc.  Equation  (2. 1-10) 
can  be  rewritten 


T 


R 


c+  6  c 
0 


(2. 1-11) 
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Since  the  uncertainty,  6c,  is  small  compared  to  cp  Eq.  (2. 1-11)  can  be 
approximated  as 


(2. 1-12) 


The  propagation  uncertainty  or  delay  due  to  uncertainties  in  the  speed  of 
light,  6T,  for  ground  waves  can  be  separated  from  the  above  equation 
as 


6T 


R6c 


(2.1-13) 


The  term  6c/pq  may  be  considered  a  random  scale  factor  error,  a. 

The  propagation  uncertainty  for  ground  waves*  can  therefore  be  modeled 
as  a  random  scale  factor,  times  the  propagation  time  of  a  wave  over  a 
similar  distance  at  the  nominal  velocity  of  light  in  the  atmosphere,  R/c^ 


(2. 1-14) 


Although  a  similar  range -dependent  effect  in  propagation  uncertainty 
occurs  with  skywaves  its  influence  is  eighed  by  random  delays  due 
to  ionosphere  altitude  variations  and  ionospheric  disturbances.  For 
skywave  propagation,  the  propagation  uncertainties  are  best  modeled 


■? - - - 

These  equations  also  apply  for  other  forms  of  direct  transmission  such 
as  satellite  radio  transmissions  passing  through  the  atmosphere. 
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as  a  random  bias  error*  independent  of  the  range  between  the  transmitter 
and  receiver. 

The  effect  of  propagation  uncertainties  for  radio  navigation 
systems  is  to  cause  errors  in  determining  the  LOP’s  on  which  the 
receiver  lies. 

2. 1. 3  Geophysical  Uncertainties 

Inertial  systems  operate  under  the  assumption  that  the  earth's 
gravity  field  (net  result  of  earth’s  gravitational  attraction  and  centripetal 
acceleration  experienced  from  the  earth's  rotation)  has  equipotential  sur¬ 
faces  which  are  elliposidal  in  shape.  **  Inertial  systems  must  compensate 
for  the  specific  force  due  to  gravity  using  knowledge  of  the  gravity  field. 
However,  many  deviations  from  the  ideal  gravity  surface  exist,  appearing 
as  variations  in  the  direction  and  magnitude  of  the  gravity  vector;  the 
former  are  referred  to  as  gravity  deflections  and  the  latter  as  gravity 
anomalies.  They  result  from  the  fact  that  the  earth  is  not  a  smooth 
homogeneous  body.  Their  effect  on  the  navigator  position  error  can 
usually  be  reduced  when  external  position  fixes  are  available. 

Gravity  anomalies  and  vertical  deflections,  are  typically 
modeled  as  first-order  markov  processes  whose  correlation  times  are 


* 

This  model  applies  when  considering  tactical  missiles  due  to  the  short 
time  of  flight.  Daily  and  seasonal  variations  do  occur  which  alter  the 
model  for  longer  time  intervals. 

**  '  V 

Even  the  harmonic  equipotential  surfaces  measured  and  used  in 
satellite  applications  do  not  account  for  the  variations  discussed  here. 
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inversely  proportional  to  vehicle  velocity.  *  Vertical  deflections  are 
assumed  to  enter  the  system  in  the  form  of  small  disturbance  inputs 
to  level  loop  accelerometers.  For  analysis  of  inertial  systems  in 
tactical  missiles,  velocity  and  position  errors  can  be  computed  on  the 
basis  of  open-loop  propagation  of  these  errors;  the  tactical  missile 
flight  time  is  short  relative  to  the  84 -min.  mode  of  inertial  navigator 
level  loops.  The  vertical  deflection  disturbances,  6,  are  assumed  to 
be  random,  zero  mean  quantities  with  autocorrelation  function: 


Because  vertical  deflections  exhibit  a  spacial  rather  than  temporal  vari¬ 
ation,  the  time  constant  is  computed  from  a  distance  correlation  Dg 
by  assuming  constant  missile  horizontal  velocity,  V: 


T  =  ~ 

V 


(2. 1-16) 


The  velocity  and  position  errors  in  each  level  loop  are  gener¬ 
ated  according  to 

ev  e6(r?)  dr?  ep  =  j  ev^  dn  (2* 1_17) 


The  analysis  here  is  a  portion  of  the  work  previously  discussed  in 
Ref.  2.  Also  considered  in  Ref.  3  was  a  second-order  markov  model 
for  the  vertical  deflecticjns  and  gravity  anomalies.  The  results  obtained 
from  the  second-order  model  were  very  similar  to  the  first-order  model, 
so  only  one  set  of  analyses  is  presented  in  this  report. 
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Mean  square  values  of  position  and  velocity  errors  are  given  by 


ev  = 
e^  = 


?.  r*3  T 

2T/4 

o  6 


[3  2 


T2te 


-t/T6 


+T' 


-t/T 


(2. 1-18) 


Under  the  assumption  of  constant  missile  velocity,  the  quantity  t/Tg  at 
the  target  is  the  ratio  of  missile  range,  R,  to  Dg.  Defining 


and  dropping  the  subscript  notation,  the  expressions  for  mean  square 
errors  at  the  target  are: 


!  ! 


I 


y 


s 

e 

v 


1 

ep 


2Ta2t[l-|(l-e-a)j;  t  =  | 

0rr  2*3  fl  1  1  -aA  1  _axl 

2Ta  X  [3  "  2a  “  2  e  3  ^  "e 


(2. 1-19) 


Figure  A.  3-1  illustrates  these  relations.  Since  external  position  fixes 
can  reset  position  errors  in  inertially  aided  systems  these  errors  are 
generally  considered  important  only  for  pure  inertial  navigation. 
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Figure  2. 1-5  Mean  Square  Normalized  Velocity  and  Position 
Errors  as  a  Function  of  the  Parameters  ”a" 

2. 1. 4  Electronic  Clock  Accuracy 

Some  radio  navigation  systems  such  as  Direct  Ranging  LORAN 
locate  the  receiver  on  a  circular  LOP  by  measuring  the  time  required  for 
a  signal  to  travel  from  a  master  station.  Such  systems  require  the  user 
to  have  a  clock  on  board  to  determine  the  time  of  signal  arrival;  informa¬ 
tion  giving  the  time  the  signal  was  sent  and  the  station  from  which  it  was 
sent  is  contained  in  the  signal.  The  propagation  time  is  therefore  the 
difference  between  the  time  of  signal  arrival  (from  the  user  clock)  and  the 
time  the  signal  was  sent.  This  can  be  converted  to  range  from  the  trans¬ 
mitting  station  by  multiplying  the  propagation  time  by  the  velocity  at  which 
radio  signals  propagate,  hi  most  cases  the  user  clock  is  a  quartz-crystal 
clock  and  this  will  be  assumed  for  the  present  studies.  The  accuracy  in 
determining  the  time  of  signal  arrival  effects  the  error  in  determining  the 
LOP  on  which  the  receiver  lies.  The  error  sources  encountered  in  crystal 
clocks  are  described  below. 
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Initial  Phase  Error  —  This  error  is  the  actual  time  difference 
between  the  receiver  clock  and  the  reference  master  station  clock  before 
any  calibration  procedures  begin.  .Since  the  phase  of  an  unsynchronized 
crystal  oscillator  is  arbitrary,  the  Initial  phase  error  may  be  quite  large 
(e.  g. ,  1-10  p  sec).  Tiie  residual  error  in  the  initial  phase  error  after 
synchronization  with  the  master  clc-k  contributes  to  the  LOP  location 
error. 


Frequency  Offset  Error  —  Frequency  offset  error  refers  to  the 
difference  in  frequency  between  the  quartz -crystal  clock  and  the  master 
clock  immediately  after  synchronization  has  been  achieved.  These  errors 
are  a  combination  of  frequency  measurement  error  and  limitations  on  the 
adjustability  of  the  quartz -crystal  clock  frequency.  Any  initial  frequency 
error  will  result  in  a  linear  phase  error  build-up  with  time.  This  phase 
error  build-up  is  determined  by  considering  the  nominal  clock  frequency 
fQ  and  tq,  the  corresponding  period.  The  frequency  fractional  error,  6, 
is  given  by 

f-f  A, 

6  =  =  (2.1-20) 
o  o 

where  f  is  the  actual  frequency.  It  can  be  shown  that  the  clock  time  error 
6T  after  t  seconds  of  operation  is 

6T  =  (2.1-21) 

o 

This  implies  that  any  fixed  frequency  error  Af/fQ  will  cause  a  system 
time  error  whose  absolute  value  is  proportional  tc  time  since  synchroni¬ 
zation. 
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Long-Term  Stability  —  Crystal  clocks  exhibit  changes  in  fre¬ 
quency  over  a  period  of  days.  These  changes  are  the  result  of  crystal 
aging  or  parametric  changes  that  are  partially  predictable.  Any  near¬ 
constant  rate  of  change  of  frequency  corresponds  to  a  second  derivative 
of  the  time  error.  This  results  in  a  quadratic  time  error  build-up. 

Short-Term  Stability  —  Short-term  stability  refers  to  the  effect 
of  small  random  clock  frequency  variations  over  intervals  that  are  short 
compared  to  the  inverse  of  the  frequency  error,  f-fQ.  .  For  tHe  short  flight 
times  oi  tactical  missiles  this  error  source  may  be  considered  a  white 
noise  frequency  error  which  is  integrated  over  the  missile  flight  time  to 
produce  a  clock  phase  error. 

Error  Due  to  Temperature  Variations  —  The  quartz -crystal 
of  the  clock  is  temperature  sensitive,  and  exhibits  a  frequency  shift  for 
a  change  in  temperature.  The  temperature  variations  can  be  assumed  to 
have  a  spatial  distribution.  The  ’’thermal  mass"  of  the  clock,  however, 
prevents  these  temperature  variations  from  affecting  clock  frequency  in¬ 
stantaneously.  The  error  due  to  temperature  variations  may  be  considered 
a  white  noise  frequency  error  for  the  short  missile  flight  times  of  interest 
and  the  temperature  variations  that  are  likely  to  be  encountered. 

Vibration -Induced  Errors  —  Quartz -crystal  clocks  exhibit  a 
shift  in  frequency  when  subjected  to  random  vibrations.  This  frequency 
shift  depends  upon  the  magnitude  and  spectrum  of  the  induced  vibration. 
Since  the  correlation  time  of  the  vibration  induced  error  is  extremely 
small,  it  may  be  shown  (Ref.  1)  that  for  typical  measurement  rates 
(Measurements  referring  to  the  measurement  of  propagation  time  from 
a  master  station)  this  error  may  be  model  as  a  white  noise  frequency 
error  which  is  integrated  to  get  clock  phase  error. 
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The  above-listed  error  sources  may  be  mathematically  repre 
sented  as  (Ref.  1). 

% 


&PC  =  %  +  “o*  +  6it2 


6<p  =  the  total  clock  phase  error  in  interval  t  seconds 
c 

&p0  =  initial  phase  error 

6f0  =  initial  frequency  offset 

=  long-term  stability  coefficient 

E(t)  =  white  noise  frequency  errors  due  to  short¬ 
term  stability,  temperature  variations  and 
vibrations. 


This  equation  may  be  modeled  as  shovn  in  Fig.  2. 1-6. 


K-4S28  a 


8,  8f0  8«p0 


E(«) 


Figure  2. 1-6  Clock  Error  Model 


In  some  modes  of  operation  the  clock  errors  may  be  calibrated  out. 
The  clock  error,  or  clock  calibration  error  if  calibration  is  used, 
directly  effect  the  error  in  the  range  measurement  for  such  systems 
as  Direct  Ranging  LORAN. 
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2.1.5  Other  External  Aid  Fix  Errors  . . 

For  some  aided  inertial  guidance  techniques  (radar  correlation, 
optical  correlation,  and  doppler),,  the  limiting  accuracy  may  be  directly 
reli  ted  to  the  fix  accuracy  of  the  aiding  device.  However,  errors  in  the 
aiding  device  occur  for  reasons  other  than  those  discussed  in  the  previous 
sections.  The  errors  in  these  aiding  techniques  are  dependent  on  the  ex¬ 
ternal  aid  and  will  be  discussed  separately  in  the  sections  which  cover 
these  guidance  techniques. 


2. 2  LIMITING  ACCURACY  OF  MISSILE  GUIDANCE  TECHNIQUES 

2. 2. 1  Pure-Inertial  Guidance 

Pure  inertial  guidance  in  tactical  missiles  refers  to  the  use  of 
an  ;  iertial  navigator  to  guide  the  missile  to  a  precomputed  target  location. 

The  accuracy  of  pure  inertial  navigation  for  use  in  tactical 
'misfiles  is  limited  by  the  errors  induced  by  vertical  deflections  and 
gravity  anomalies.  The  rms  magnitude  of  these  errors  may  be  com¬ 
puted  using  Fig.  2. 1-5. 

To.  illustrate  the  use  of  the  figure  the  following  example  is 

-3 

©resented.  Representative  values  for  D.  and  o,  are  30  nm  and  1.  3  x  10 
ft/uee  (8  sec  deflection)  respectively.  When 

R  =  90  nm 
V  =  10^  ft/sec 
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Figure  2. 1-5  indicates  that 

q  ss  0. 47  ft/ sec 
ev 

oe  a-  138  ft 

per  level  loop,  at  the  target. 

hi  the  vertical  direction,  gravity  anomalies  (variations  in  the 
magnitude  of  gravity)  cause  navigation  errors.  These  random  variations 
can  be  modeled  in  a  manner  analogous  to  deflections  of  the  vertical.  The 
vertical  position  loop  of  an  unaided  inertial  navigator  is  unstable.  How¬ 
ever,  the  missile  flight  time  is  so  short  that  the  instability  hardly  has  time 
to  contribute  error  growth.  Considering_only  the  open  loop  vertical  navi¬ 
gation  errors,  the  equations  for  e$  and  e^  are  the  same  as  those  provided 
above.  Only  the  magnitude  of  the  correlation  time  T  and  disturbance 
standard  deviation  a  may  differ  from  those  of  vertical  deflections. 
Representative  values  of  D  and  or  for  gravity  anomalies  are  30  nm  and 
3  x  10“5  g’s.  The  resulting  rms  velocity  and  position  errors  in  the  ver¬ 
tical  direction  are  about  0.35  ft/sec  and  103  ft  at  the  end  of  a  90-nm 
flight  (V  =  lO^ft/sec),  Thus,  the  rss  velocity  and  position  errors  due  to 
vertical  deflections  and  gravity  anomalies  are  0.75  ft/sec  and  220  ft  in 
this  example.  Errors  of  this  magnitude  are  not  serious  when  a  target 
homing  device  is  used  for  terminal  guidance  systems  but  serve  to  illus¬ 
trate  an  accuracy  limit  for  pure  inertial  missile  guidance  systems  if  no 
external  position  and  velocity  aids  are  available.  Similar  calculations 
can  be  carried  out  for  other  missile  ranges  or  velocities. 

2. 2. 2  LORAN/lnertial  Guidance 

Conventional  LORAN  determines  the  LOP’s  on  which  the  re¬ 
ceiver  lies  by  measuirng  the  difference  in  the  time  of  arrival  of  signals 
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from  two  stations.  The  LOP’s  are  thus  hyperbolic  in  shape.  It  is  as¬ 
sumed  in  this  analysis  that  the  receiver  location  measurements  provided 
by  LORAN  are  optimally  combined  to  estimate  position  errors  in  the 
inertial  navigator.  If  many  LORAN  fixes  are  taken  then  all  white  noise 
error  sources  contributing  to  the  LORAN  fix  accuracy  (such  as  receiver 
noise)  are  averaged  out.  Random  bias  errors  in  the  LORAN  measure¬ 
ments  will  determine  the  ultimate  limits  on  accuracy.  These  bias  errors 
cannot  be  estimated  since  the  distance  of  the  missile  from  the  LORAN 
stations  is  typically  large  compared  to  the  range  of  the  missile;  little 
geometry  change  occurs  and  the  bias  type  errors  in  the  measurements 
are  not  observable. 

The  two  m~jor  factors  whic'i  limit  the  accuracy  of  a  conventional 
LORAN  system  are  propagation  uncertainties,  which  cause  errors  in  de¬ 
termining  the  LOP’s,  and  geometric  dilution  which  magnifies  radial  errors 
in  situations  of  poor  geometry.  To  illustrate  how  the  equations  derived 
in  Section  2. 1. 1  may  be  used  to  determine  the  accuracy  of  conventional 
LORAN  an  example  is  presented  below. 

Assume  that  a  missile  is  flying  in  a  configuration  of  LORAN 
stations  as  shown  in  Fig.  2.  2-1.  As  stated  previously  due  to  the  short 
flight  range  of  the  missile  compared  to  the  range  to  the  stations,  the 
geometry  may  be  assumed  to  be  fixed. 

Conventional  LORAN  measures  the  difference  in  time  of  ar¬ 
rival  of  signals  from  two  stations.  LORAN  waves  are  ground  waves* 

* 

Skywave  propagation  becomes  important  only  at  large  distances  from 
LORAN  stations.  At  these  ranges  the  accuracy  greatly  decreases  so 
that  navigation  at  large  distances  from  LORAN  stations  will  not  be 
considered. 
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l 

L 

L 


The  error  in  the  time  difference  measurement,  AT  may  be 
viewed  as  the  range  from  one  station  being  in  error  by  an  amount  cQ  AT/2 
and  the  range  to  the  other  station  by  an  amount  -c0AT/2.  Since  these 
range  errors  are  small  compared  to  the  range  to  each  station,  the  indi¬ 
cated  LOP  can  be  considered  shifted  by  an  amount  c  as  shown  in 
Fig.  2. 2-2.  From  the  geometry,  the  error  in  determining  the  LOP  is 


c  = 


c  AT 
o 


2sin(<p/2) 


(2.  2-2) 


where  <p  is  the  angle  between  the  vectors  from  the  leceiver  to  the  sta¬ 
tions  (see  Fig.  2. 2-2).  Again  using  subscripts  to  denote  stations, 
Eqs.  (2. 2-2)  and  (2. 2-3)  can  be  combined  to  give 


1,2  2  sin(<pj/2) 


B20!2"R3a3 
F2,3  =  2sin(«P2/2) 


(2. 2-3) 


Squaring  these  equations  and  taking  the  expected  value  gives 


7~\  __l 
1,2  , 
4  ebv 


^  («i  <4 + R?  2r2ri  °i) 


7~. 

2,3 


4  sin 


\  (R2  a2  +  R3  “3  “  2R2R3  a2aV 

V2  / 


! 

I 
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Figure  2.  2-2  LOP  Error  Generated  from  Error  in  Propagation 
Time  Difference  Measurement  from  Two  Stations 
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Assuming  the  LORAN  propagation  uncertainties  from  different  stations 
are  uncorrelated,  Eq.  (2. 2-4)  can  be  reduced  to 


(2. 2-5) 


where 


ctj  =  standard  deviation  in  the  LOP  position  error 
1  from  time  difference  measurement  from 
stations  1  and  2 

=  standard  deviation  in  the  LOP  position  error 
from  time  difference  measurement  from 
stations  2  and  3 

p  «  correlation  coefficient  between  the  position 
errors  for  the  two  LOP’s 

i 

cr  ,o  ,o  =  the  standard  deviations  of  the  propagation 
1  “3  scale  factor  error  for  each  of  the  three 

stations  .. 
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The  standard  deviation  of  the  radial  receiver  location  error, 
od>  can  be  found  by  combining  Eqs.  (2.  2-5)  and  (2. 1-8).  To  illustrate 
the  use  of  these  equations  in  solving  for  the  limiting  accuracy  (radial 
error)  for  the  LQRAN/inertial  system  the  following  example  is  pre¬ 
sented:  The  station  geometry  shown  in  Fig.  2.2-3  was  assumed.  The 
receiver  location  was  varied  about  this  station  configuration  and  Eq. 

(2. 2-5)  was  evaluated  to  determine  the  standard  deviations  of  the  LOP 
position  errors  and  the  correlation  coefficient  between  them.  The  stan¬ 
dard  deviation  of  the  radial  error  was  then  found  from  Eq.  (2. 1-8). 
These  equations  were  evaluated  assuming  the  propagation  uncertainty 
scale  factors,  o^,  have  a  magnitude  of  0.  0001.  The  results  are  pre¬ 
sented  in  Fig.  2.  2-4.  This  plot  provides  contour  lines  of  equal  radial 
error.  It  is  interpreted  as  meaning  that  the  limiting  receiver  location 
error  (radial)  is  the  value  given  for  the  contour  line  when  the  receiver 
lies  on  that  line.  Interpolation  may  be  used  to  determine  radial  location 
errors  when  the  receiver  lies  between  contour  lines. 

R-78S9 

STATION  1 


STATION  3 

Figure  2.  2-3  LORAN  Station  Geometry  Assumed 
in  the  Example 
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O  TRANSMITTER 

PROPAGATION  ANOMALY  SCALE  FACTOR  MO*4  , 


RADIAL  ERROR 
CONTOUR  LEVEL  VALUES 


A  = 

1000  ft 

F  = 

10,000  ft 

B  * 

1000  ft 

G  « 

io,ooo  ft 

C  ? 

3000  ft 

H* 

30,000  ft 

D  * 

5000  ft 

1  a 

40,000  ft 

E  * 

7500  ft 

J  * 

50,000  ft 

Figure  2.  2-4  Limiting  Radial  Error  for  LORAN/lnertial 
Navigation:  10 "4  Propagation  Anomaly 
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The  magnitude  of  the  propagation  uncertainty  scale  factor,  o  , 
is  dependent  on  the  knowledge  of  the  local  propagation  variations.  It  has 
been  established  that  these  may  be  determined  in  a  given  area  to  the  extent 
that  the  residual  error  is  given  by  =  0. 0001.  However  without  pre¬ 
surveying,’  aa  may  run  as  high  as  0. 01.  Examination  of  the  equations  used 
to  obtain  Fig.  2. 2-4  shows  that  the  results  may  be  linearly  scaled  to  ac¬ 
count  for  other  values  of  the  propagation  uncertainty  scale  factor.  For 
example,  if  the  propagation  scale  factor  is  0. 01  the  values  of  the  contour 
levels  are  simply  multiplied  by  100. 

To  illustrate  a  particular  example  assume  that  the  missile  is 
located  700  nm  from  Station  2  and  equidistant  from  Stations  1  and  3 
(approximately  600  nm  from  each).  From  the  figure  it  may  be  shown 
that  the  limiting  accuracy  for  the  LORAN/Inertial  system  at  this  location 
is  approximately  2800  ft  for  a  propagation  scale  factor  of  0. 0001.  This 
would  be  acceptable  for  missile  navigation  if  a  homing  sensor  is  used  for 
final  guidance  to  the  target.  However  if  there  is  no  knowledge  of  propaga¬ 
tion  anomalies  and  o  =  0. 01,  the  missile  location  error  is  280, 000  ft, 
which  is  probably  unacceptable  even  with  a  homing  phase.  Similar  calcula¬ 
tions  may  be  carried  out  to  find  the  radial  error  for  other  missile  station 
configurations. 


With  Direct  Ranging  LORAN  (DRL)  the  receiver  location  is 
determined  on  a  circular  LOP  by  measuring  the  time  required  for  a  sig¬ 
nal  to  travel  from  a  master  station.  This  is  accomplished  by  comparing 
the  time  the  LORAN  signal  was  sent  (information  giving  the  time  the  sig¬ 
nal  was  sent  and  the  station  from  which  it  was  sent  is  contained  in  the 
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signal)  and  the  time  of  signal  arrival  (determined  from  an  onboard  user 
clock).  This  time  difference  can  be  converted  to  range  from  the  trans¬ 
mitting  station  by  multiplying  the  propagation  time  by  the  velocity  at 
which  radio  signals  propagate.  Measurements  from  only  two  stations 
are  required  to  locate  the  receiver  (assuming  an  independent  measure 
of  altitude),  but  measurements  from  three  stations  are  typically  used 
with  the  redundant  information  being  used  to  calibrate  errors  in  the  user 
clock  and  estimate  propagation  delays.  It  is  assumed  in  this  analysis 
that  the  receiver  location  measurements  from  three  stations  are  optimally 
combined  to  estimate  position  errors  in  the  inertial  navigator,  calibrate 
the  user  clock,  and  estimate  propagation  delays.  If  many  LORAN  fixes 
are  taken  then  all  white  noise  error  sources  contributing  to  the  LORAN 
fix  accuracy  (such  as  receiver  noise)  are  averaged  out.  The  measurement 
geometry  (the  measurement  matrix  in  the  Kalman  filter  formulation)  may 
be  assumed  fixed  since  the  distance  of  the  missile  from  the  LORAN  sta¬ 
tions  is  typically  large  compared  to  the  range  of  the  missile;  little  geo¬ 
metry  change  occurs.  It  is  also  assumed  that  the  navigator  position  errors 
are  very  large  before  the  beginning  of  the  processing  of  LORAN  fixes 
thus  making  the  limiting  accuracy  of  the  DRL/inertial  system  dependent 
on  only  the  LORAN  fix  accuracy  and  not  the  initial  navigator  position 
errors.  Since  the  time  of  flight  of  a  tactical  missile  is  short  the  time- 
■  varying  errors  in  the  inertial  system  and  user  clock  are  ignored. 

Since  the  navigator  error  dynamics  are  not  considered,  the 
important  errors  which  influence  the  limiting  accuracy  of  the  system  are:* 

* 

Notice  that  this  list,  together  with  assumption  stated  immediately  before 
it,  precludes  any  inertial  system  error  source  from  having  an  effect  on 
the  limiting  accuracy  of  the  DRL/inertial  system.  Tne  key  here  is  that 
we  are  computing  limiting  accuracy  rather  than  practical  accuracy. 
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Xj  =  a  navigator  position  error* 

Xg  =  b  navigator  position  error* 
x,j  =  user  clock  error 

'x^  =  propagation  anomaly  scale  factor  from  station  1 

Xg  =  propagation  anomaly  scale  factor  from  station  2 

Xg  =  propagation  anomaly  scale  factor  from  station  3 

We  designate  x  as  the  state  vector  of  these  errors 

*T  =[xl  *2  X3  X4  *5  \]  (2'2-C) 

and  consider  Pq  to  be  the  initial  covariance  matrix  of  these  errors  before 
any  updates  using  the  LORAN  information.  The  covariance  matrix  Pn 
after  n  sets  of  updates  is  given  by  Eq.  (2.  2-7).' 

-1  T  -1 
P*  +  H1R  lE 
o 

-1  T  -1 
PQ  +  2ITR  H 

P“14nHTR"1H  (2.2-7) 

The  a,b  coordinate  f"ame  is  any  orthogonal  set  of  horizontal  coordi¬ 
nates. 
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where  R  is  the  measurement  noise  matrix  for  the  LORAN  fh  '  s  and  II  if. 
the  measurement  matrix.  *  By  appropriate  matrix  manipulr  J  ,n  fcq. 

(2. 2-7)  may  be  written 

Pa  =  PQ  -  P0HT(HP0HT+R/n)  HPq  (2.2-0) 

As  the  number  of  fixes  approaches  infinity  (the  limiting  case)  the  equation 


reduces  to 


p  =  P  -PHJ 

*00  o  o 


KhT)" 


(2.2-9) 


The  covariance  of  the  radial  error  can  then  be  easily  found 


from  the  elements  of  as 


cov  (radial  error)  =  cov  (r^)  +  cov  (xg)  (2. 2-10) 

where  and  Xg  are  defined  in  the  state  vector. 

For  a  general  three -station  configuration  as  shown  in  Fig.  2.2-5 
the  measurement  matrix,  H,  can  be  defined  for  the  state  vector  in  Eq. 

(2.  2-6)  as  follows: 

'(Val)  (VM  1  E.  0  o’ 

RX  Ri  1 


Rx  0  0 


(ar"a2)  (Vb2)  1 

—  _ ""  1  n  ^ 


0  r2  ° 


(2.  2-11) 


(Va3)  (br~b3)  t 


0  0  R< 


*This  is  constant  due  to  the  assumption  of  invariant  fix  geometry. 
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To  illustrate  the  use  of  the  above  method  for  finding  the  limit¬ 
ing  accuracy  (radial  error)  for  the  DRL/inertial  system  the  following 
example  is  presented:  The  statio.n  geometry  shown  in  Fig.  2. 2-3  was 
assumed.  The  receiver  location  was  varied  about  the  station  configura¬ 
tion  and  Eq.  (2. 2-11)  was  evaluated  to  determine  the  measurement  mat¬ 
rix,  H,  at  each  receiver  location.  Equation  (2.  2-9)  was  solved  for  each 
location  to  determine  the  covariance  of  the  location  errors  after  an  infinite 
number  of  fixes.  *  The  covariance  of  the  radial  error  may  then  be  found 
from  Eq.  (2. 2-10).  Several  choices  of  the  initial  covariance  matrix,  P  , 
were  considered  to  demonstrate  the  radial  error  sensitivity  to  variations 
in  values  of  user  clock  accuracy  and  propagation  anomalies.  Table  2. 2-1 
summarizes  the  rms  values  used  in  the  initial  covariance  matrix  for  this 
example  and  provides  the  figure  number  in  which  the  results  are  plotted. 
These  plots  provide  contour  lines  of  equal  radial  error.  They  may  be 
interpreted  in  the  same  manner  as  the  one  for  the  LORAN/inertial  system 
and  interpolation  may  be  used  between  curves. 

As  stated  previously,  the  initial  navigator  position  errors  are 
assumed  to  be  very  large  to  eliminate  the  dependence  of  navigator  accu¬ 
racy  from  the  problem.  This  is  representative  of  the  case  in  which  the 
navigator  has  poor  prior  knowledge  of  the  target  location;  i.  e. ,  only 
the  LORAN  grid  coordinates  of  the  target  are  known  and  the  receiver 
must  be  located  in  the  LORAN  grid. 


As  may  be  seen  from  Eq.  (2.  2  -8)  the  actual  number  of  fixes  needed  to 
approach  this  steady -state  solution  is  dependent  on  the  magnitude  of  the 
measurement  noise  matrix,  R.  For  typical  systems  the  number  of  fixes 
needed  to  reach  steady  state  is  small  and  easily  achieved  under  operational 
conditions. 
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TABLE  2.  2-1 

SUMMARY  OF  INITIAL  COVARIANCE  VALUES 
FOR  EXAMPLE  RUNS 


f~  Initial  RMS  Value 

Figure  Number 

Showing  Limiting  Accuracy 
(radial  error) 

User  Clock 
Error  (nsec) 

Propagation  Anomaly 
Scale  Factor  Error* 

300 

10“4 

2.2-6 

500 

10-4 

2.2-7 

300-500 

10"2 

2.2-8 

User  clock  errors  of  300  nsec  are  considered  typical  values 
for  quartz-crystal  clocks  while  500  nsec  is  considered  a  worst  case 
(Ref.  1).  Therefore,  both  situations  were  studied.  Also  various  values 
for  the  propagation  anomaly  scale  factor  were  used  since  the  magnitude 
is  dependent  on  the  knowledge  of  the  local  propagation  variations.  It  has 

bean  established  (Ref.  1)  that  these  may  be  de  termined  in  a  given  area 

-4 

to  the  extent  that  the  residual  scale  factor  error  is  10  .  However, 

-2 

without  presurveying  values  may  run  as  high  i.a  10  .  Both  cases  were 
considered. 

Comparison  of  Figs.  2.  2-6  and  2.  2-7  shows  that  the  user  clock 
error  has  its  most  significant  influence  on  the  radial  error  sensitivity 
near  the  LORAN  transmitters  when  propagation  anomalies  are  small. 

This  occurs  because  when  the  missile  is  close  to  the  transmitters  the 
effects  of  propagation  uncertainties  are  small  due  to  the  short  distances 


* 

The  propagation  delays  to  the  receiver  from 
assumed  uncorrelated  (Ref.  1). 


different  stations  are 
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the  radio  waves  must  travel.  At  large  distances  the  propagation  delays 
become  the'  largest  contributor  to  the  radial  error,  thus  decreasing  the 
sensitivity  to  user  clock  error. 

For  large  magnitudes  of  the  propagation  scale  factor  error 
(areas  which  are  not  presurveyed)  the  propagation  uncertainties  are  the 
dominant  error  sources  contributing  to  the  Direct  Ranging  LORAN/inertial 
system.  Therefore  Fig.  2.2-8  applies  for  the  entire  range  of  user  clock 
errors  from  300  to  500  nsec.  When  the  receiver  is  located  near  the  cen¬ 
ter  of  this  figure,  within  the  contour  line  of  lowest  radial  error  (50,000  ft), 
the  location  accuracy  is  approximately  50, 000  ft.  The  accuracy  does  not 
increase  as  the  receiver  moves  nearer  the  transmitting  stations  since  the 
fix  geometry  (crossing  angles  of  the  LOP's)  becomes  less  favorable  in  this 
area  for  estimation  of  the  propagation  delays. 

It  is  interesting  to  compare  the  accuracy  of  the  Direct  Ranging 
LORAN/inertial  systems  with  that  of  hyperbolic  LORAN/inertial  systems, 
reported  in  the  previous  section,  hi  that  section  the  example  chosen  as¬ 
sumed  the  same  station  configuration  as  above  with  the  receiver  located 
700  nm  from  Station  2  Eind  equidistant  from  Stations  1  and  3  (approximately 
600  nm  from  each).  The  results  showed  that  with  hyperbolic  LORAN/ 
inertial  the  limiting  accuracy  at  this  point  was  2, 800  ft  in  presurveyed 
areas  (propagation  scale  factor  =  10"^).  From  Fig.  2.  2-7  it  is  found  that 
for  the  worst  case  user  clock  accuracy  the  DRL/inertial  radial  error  is 
approximately  760  ft.  Thus  an  improvement  of  at  least  a  factor  of  3.  6  is 
achieved  with  Direct  Ranging  LORAN.  This  results  from  better  fix  geo¬ 
metry  with  Direct  Ranging  LORAN  (circular  LOP’s  instead  of  hyperbolic) 
and  the  fact  that  three  measurements  are  made  which  are  redundant  and 
therefore  user  clock  accuracy  and  propagation  uncertainties  may  be  esti¬ 
mated  to  a  certain  extent. 
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Examination  of  Jigs.  2. 2-6,  2. 2-7,  and  2. 2-8  shows  that  the 
DRL/inertial  technique  provides  adequate  accuracy  for  midcourse  guid¬ 
ance  in  presurveyed  areas  but  probably  is  not  acceptable  for  terminal 
guidance  without  additional  terminal  homing  devices.  In  unsurveyed  areas 
the  technique  is  generally  not  acceptable. 

2. 2. 4  OMEGA/lnertial  Guidance 

OMEGA  is  a  long  range  radio  navigation  system  whose  accuracy 
is  limited  by  propagation  uncertainties,  OMEGA  measures  the  difference 
in  phase  of  signals  from  two  stations  to  locate  the  receiver  on  a  hyperbolic 
line  of  position.  The  location  of  the  receiver  is  then  the  intersection  of 
two  LOP's.  As  with  LORAN,  the  receiver  location  error  is  magnified 
for  situations  of  poor  geometry  (geometric  dilution). 

We  will  assume  that  the  receiver  location  measurements  pro¬ 
vided  by  OMEGA  are  optimally  combined  to  estimate  position  errors  in 
the  inertial  navigator.  As  with  LORAN  little  geometry  change  occurs 
since  the  range  of  a  tactical  missile  is  typically  small  compared  to  the 
distance  of  the  missile  from  the  stations.  For  similar  reasons  as  those 
presented  in  the  previous  section  it  can  be  shown  that  the  accuracy  of  the 
OMEGA/inertial  system  can  be  related  to  the  random  bias  errors  in  the 
OMEGA  measurements  (all  white  noise  error  sources  may  be  averaged 
out  by  taking  many  measurements). 

OMEGA  signals  are  the  skywave  type  (see  Section  2. 1. 2)  and 
the  error  in  determining  a  line  of  position  is  considered  to  be  a  random 
bias  which  is  undorrelated  with  bias  errors  in  LOP  determination  from 
"other  combination*,  of  stations.  The  present  accuracy  of  OMEGA  provides 
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an  error  in  determining  LOP's  with  a  standard  deviation  of  from  0. 75  to 
l,0nm  (Ref.  4).  It  is  estimated  that  the  theoretical  limit  on  this  accu¬ 
racy  is  0. 5  nm. 

Since  the  precision  of  OMEGA  is  geometry  dependent  an  example 
will  be  given  to  illustrate  the  calculations  necessary  to  determine  the  limit¬ 
ing  accuracy.  For  this  purpose  of  this  example  an  OMEGA  station  configu¬ 
ration  as  illustrated  in  Fig.  2.  2-9  is  assumed.  This  is  representative 
of  the  three  stations  located  in  Aldra,  Norway;  Rome,  New  York;  and 
Trinidad,  West  Indies.  The  LOP’s  are  determined  by  the  phase  difference 
of  signals  from  stations  1  and  2  and  from  stations  2  and  3.  Since  the 
measurement  errors  are  uncorrelated-and  the  standard  deviation  of  their 
magnitudes  should  be  approximately  equal,  Eq;  (2. 1-8)  can  be  simplified 
to  give  a  standard  deviation  for  the  radial  error  as  follows: 


(2.2-13) 


where  cr  is  the  standard  deviation  of  the  bias  errors  in  the  measure  - 
m 

ment  of  LOP’s  and  and  <p2  are  the  angles  shown  in  Fig.  2. 2-9.  The 
receiver  location  was  varied  about  the  station  configuration  and  Eq. 

(2. 2-13)  was  evaluated.  The  limiting  accuracy  (radial  error)  for  OMEGA/ 
inertial  navigation  determined  in  this  manner  is  presented  in  Fig.  2. 2-10 
for  a  1  nm  propagation  anomaly  bias  error.  This  figure  presents  con¬ 
tour  lines  of  equal  radial  error  and  may  be  interpreted  in  the  same  man¬ 
ner  ae  the  one  for  the  LORAN/inertial  system.  Interpolation  may  be 
used  between  the  contour  lines.  Examination  of  Eq.  (2. 2-13)  shows  that 
the  contour  level  values  may  also  be  linearly  scaled  to  evaluate  other 
values  of  the  propagation  anomaly  bias  error.  For  example  if  a 
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STATION  2  STATION  3 

Figure  2.2-9  Configuration  of  Missile  Flight  Geometry 

propagation  anomaly  bias  error  is  0.  5  nm  the  contour  level  values  are 
multiplied  by  one -half. 

As  may  be  seen  from  Fig.  2. 2-10,  over  a  range  of  propagation 
anomaly  bias  errors  from  0.  5  to  1.0  nm  the  accuracy  of  OMEGA/inertial 
navigation  is  probably  acceptable  for  midcourse  guidance  in  areas  of 
good  station/receiver  geometry.  However,  a  homing  sensor  for  termi¬ 
nal  guidance  would  be  needed.  Although  these  results  apply  to  the  parti¬ 
cular  station  configuration  shown  in  Fig.  2. 2-9,  similar  calculations 
can  be  carried  out  for  other  station  configurations. 

2. 2. 5  DME/lnertial  Guidance 

Distance  Measuring  Equipment  (DME)  systems  are  typically 
short  range  radio  navigation  systems  which  measure  range  or  range  and 
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Figure  2.2-10  Limiting-Radial  Error  for  OMEGA/lnertial 
Navigation;  1  nm  Propagation  Anomaly 
Bias  Error 
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velocity  of  the  receiver  from  one  or  more  transponders.  The  accuracy 
of  DME  systems  is  a  function  of  the  number  of  transponders,  flight  path, 
and  sophistication  of  the  data  processing  techniques  used*  (number  of 
error  sources  modeled).  Due  to  the  unlimited  number  of  flight  paths 
and  station  configurations  possible,  only  a  qualitative  discussion  of  the 
accuracy  of  some  particular  configurations  will  be  given.  Quantitative 
results  would  require  a  separate  analysis  for  each  flight  configuration. 

One  Transponder  Configuration  —  A  typical  one  transponder 
configuration  is  shown  in  Fig.  2. 2-11.  With  this  configuration,  if  the 
relative  position  between  the  transponder  and  missile  is  approximately 
known  initially  by  the  inertial  navigator,  as  the  geometry  changes  along 
the  flight  path  the  position  and  velocity  of  the  missile  with  respect  to  the 
transponder  can  be  estimated  in  the  two  horizontal  dimensions  (x  and  y). 
To  perform  this  location  in  the  two  horizontal  dimensions  an  independent 
knowledge  of  altitude  are  required.  For  this  configuration  the  horizontal 
errors  in  missile,  location  are  highly  dependent  on  the  altitude  error. 

By  properly  taking  into  account  bias  or  scale  factor  errors 
due  to  propagation  uncertainties  or  other  measurement  errors  in  the 
Kalman  type  filter  used  for  data  processing  it  is  possible  to  estimate 
and  remove  their  effects.  This  occurs  since  the  effect  of  bias  and 
scale  factor  errors  change  with  the  varying  geometry  along  the  flight 
path. 


As  stated  above  a  necessary  condition  for  estimating  position 
and  velocity  is  some  knowledge  of  relative  position  between  the  trans¬ 
ponder  and  navigator;  the  navigator  position  errors  become  important 

.  •y— —  '  * 

For  the  discussions  which  follow  optimal  data  processing  is  assumed 
with  measurements  being  incorporated  by  a  Kalman  type  filter. 
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Figure  2.2-11  One -Transponder  Flight  Configuration 

in  the  ability  to  resolve  the  range  measurement  into  two  components, 
x  and  y.  Therefore  the  navigator  position  errors  enter  into  the  equations 
used  for  estimating  these  very  quantities.  Typically  no  difficulty  is  en¬ 
countered  as  long  as  the  position  errors  are  a  small  fraction  of  the  in¬ 
dicated  range  to  the  target.  However  it  is  necessary  to  simulate  the 
system  to  verify  this  fact  for  a  particular  flight  path  and  initial  set  of 
navigator  position  errors,  since  the  dependence  of  the  measurements 
on  the  position  errors  invalidates  normal  covariance  type  analyses  (the 
position  errors  no  longer  maintain  their  zero  mean  property). 

In  summary  the  accuracy  of  a  one  transponder  DME/inertial 
system  may  be  said  to  depend  on  the  following  factors: 

•  Knowledge  of  the  missile  altitude. 

•  Initial  knowledge  of  the  relative  position 
between  transponder  and  missile. 
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•  The  sophistication  of  the  filter  in  estimating 
bias  or  scale  factor  errors  due  to  propagation 
uncertainty  and  other  range  or  velocity  errors. 

•  Favorability  of  the  flight  geometiy  with  respect 
to  the  transponder. 

•  Two  Transponder  Configurations  —  A  representative  two 
transponder  configuration  is  shown  in  Fig.  2.  2-12.  Each  transponder 
measures  a  range  to  the  receiver  locating  it  on  a  sphere  about  the  trans¬ 
ponder.  The  intersection  of  the  two  spheres  positions  the  receiver  on  a 
circle.  A  knowledge  of  altitude  is  necessary  to  determine  the  receiver 
location  on  this  circle..  Again  the  altitude  error  is  a  significant  factor 
in  determining  the  horizontal  location  accuracy.  The  two  transponder 
configuration  does  eliminate  the  need  for  a  priori  coarse  knowledge  of 
relative  missile  and  transponder  positions. 

As  with  the  single  station  configuration,  by  properly  modeling 
bias  or  scale  factor  errors  in  the  transponder  measurements  their  effect 
may  be  estimated  and  removed  due  to  geometry  changes  along  the  flight 
path.  This  would  add  to  the  data  processing  complexity  but  may  allow  a 
significant  increase  in  location  accuracy. 

The  two  transponder  configuration  has  an  advantage  over  the 
one  transponder  system  discussed  above  in  that  the  additional  range 
measurement  (from  the  second  transponder)  will  both  increase  the  sys¬ 
tem  accuracy  and  decrease  the  estimation  time  of  vital  parameters 
(position,  velocity,  etc. )  for  most  flight  configurations.  For  a  further 
discussion  on  the  evaluation  of  one  particular  two  transponder  DME  sys¬ 
tem  refer  to  Ref.  3. 

Three  Transponder  Configuration  —  A  good  three  transponder 
configuration  for  missile  location  purposes  is  shown  in  Fig.  2.  2-13. 
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Since  the  missile  may  be  located  by  each  transponder  on  a  sphere  about 
itself,  the  intersection  of  the  three  spheres  defines  the  receiver  location 
in  three  dimensions.  Thus  receiver  location  does  not  depend  on  external 
knowledge  of  missile  altitude  or  prior  knowledge  of  relative  location  of 
the  navigator  and  transponders.  A  flight  path  passing  directly  over  a 
transponder  as  in  Fig.  2. 2-13  is  one  of  the  most  desirable  for  accurate 
altitude  determination.  This  configuration  allows  for  a  measure  of  alti¬ 
tude  as  the  missile  passes  directly  over  the  transponder.  There  may 
be  many  three  transponder  configurations  which  are  unfavorable  to 
altitude  determination. 

As  with  the  one  or  two  transponder  configurations  changing 
geometries  along  the  flight  path  make  it  possible  to  estimate  bias  or 
scale  factor  errors  in  the  measurement  if  these  are  properly  modeled 
in  the  filter. 
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Figure  2.  2-13  Three-Transponder  Flight  Configuration 


The  addition  of  the  third  transponder  can  provide  improved 
system  accuracies  and  decreased  estimation  times  of  vital  system  para¬ 
meters.  The  addition  of  more  transponders  over  three  provides  redundant 
measurements  and  thus  also  serves  to  improve  system  accuracies. 

2. 2. 6  Satellite/lnertial  Guidance  ' 

A  satellite/inertial  system  (for  example,  the  synchronous 
62 IB  satellite  system  in  Refs.  5  and  6)  which  determines  the  receiver 
location  on  spheres  about  each  satellite  by  measuring  the  time  required 
for  a  signal  to  travel  from  the  satellite  to  the  receiver  is  considered  for 
this  analysis.  This  is  accomplished  in  a  manner  similar  to  the  Direct 
Ranging  LORAN  system  discussed  previously  by  comparing  the  indicated 
time  on  a  satellite  clock  at  which  the  signal  was  sent  (information  giving 
the  time  the  signal  was  sent  and  the  satellite  making  the  transmission  is 


2-44 


THE  ANALYTIC  SCIENCES  CORPORATION 


contained  in  the  signal)  and  the  time  of  signal  arrival  (determined  from 
an  onboard  user  clock).  This  time  difference  can  be  converted  to  range 
from  the  transmitting  satellite  by  multiplying  the  propagation  time  by 
the  velocity  at  which  radio  signals  propagate.  Errors  are  introduced  into 
the  range  measurements  due  to  an  uncertainty  in  the  knowledge  of  the 
velocity  of  propagation. of  the  radio  waves.  Propagation  velocity  variations 
result  from  uncertainties  in  atmospheric  conductivity  caused  by  atmo¬ 
spheric  disturbances,  water  vapor  concentration,  etc.  Measurements 
from  three  satellites  are  required  to  locate  the  receiver  in  three  dimen¬ 
sions  but  a  fourth  measurement  is  typically  used  to  calibrate  errors  in 
the  satellite  and  user  clocks  and  estimate  radio  wave  propagation  uncer¬ 
tainties.  Further  details  and  results  of  the  limiting  accuracy  studies  for 
the  satellite/inertial  guidance  techniques  are  presented  in  Volume  n 
(classified  SECRET)  of  this  report,  Ref.  7. 


2.2.7  Radar  Correlation/lnertial  Guidance 


With' radar  correlation  techniques  a  position  fix  is  made  by 
correlating  a  predetermined  reference  map  or  trace  (may  be  made  from 
reconnaissance  or  satellite  photos)  with  a  "live”  return  map  generated 
by  the  missile  radar  as  it  flights  over  the  terrain.  The  location  of  the 
reference  map  is  known  and  when  the  radar  return  map  matches  this 
reference  map  the  missile  position  with  respect  to  the  reference  map  is 
known.  Further  details  and  results  of  the  limiting  accuracy  studies  for 
radar  correlation/inertial  guidance  techniques  are  presented  in  Volume  II 
(classified  SECRET)  of  this  report,  Ref.  7. 
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2. 2. 8  Optical  Corr elation/inertial  Guidance 

Optical  correlation  is  similar  in  principle  to  radar  correlation 
discussed  in  the  previous  section  except  that  optical  pictures  are  matched 
instead  of  radar  maps.  One  optical  correlation  device  is  the  NAFI  (Naval 
Avionics  Facility,  Indianapolis)  Scene  Matching  Area  Correlator  (SMAC) 
which  is  a  passive  optical  correlation  device  developed  to  provide  position 
fixes  for  a  tactical  air -to -surface  missile.  The  device  correlates  a  "live” 
image  of  the  overflown  terrain  with  a  stored  image  obtained  by  photo  re¬ 
connaissance.  Operation  of  SMAC  is  illustrated  by  Fig.  2.  2-14.  A  refer¬ 
ence  image  of  a  narrow  rectangular  strip  of  land  is  prepared  from  aerial 
reconnaissance  photographs.  This  strip  may  be  as  narrow  as  2,000  ft. 
along  the  direction  of  flight  but  its  cross -track  dimension  must  be  large 
enough  that  the  missile  will  be  almost  certain  to  fly  over  some  point  on 
the  strip.  A  negative  of  the  reference  image  is  mounted  on  a  rotating 
drum,  the  axis  of  which  is  usually  parallel  to  the  nominal  flight  path. 

A  live  image  of  the  overflown  terrain  is  focused  on  the  moving  transpar¬ 
ency.  As  the  missile  crosses  over  the  strip  of  land  photographed  the 
pre-stored  reference  image  is  scanned  over  the  live  image  by  rotating 
the  drum.  Since  the  reference  is  a  negative  image  and  the  live  scene 
represents  a  positive  image,  there  will  be  a  sharp  decrease  in  the 
amount  of  light  passing  through  the  reference  when  the  two  images  coin¬ 
cide.  This  decrease  is  due  to  a  cancellation  of  the  clear  and  opaque  areas 
:  on  the  stored  image  with  corresponding  dark  and  light  area  of  the  live 
image.  Crosstrack  position  is  determined  by  measuring  the  distance 
from  one  end  of  the  photographic  image  to  the  point  where  match  occurs. 
This  is  done  with  the  aid  of  a  magnetic  pickoff  which  senses  the  end- 
.  point  of  the  filmed  image  or.  the  rotating  drum.  The  along  track  posi¬ 
tion  is  established  by  the  fact  that  the  best  match  occurs  when  the  mis¬ 
sile  is  over  the  center  of  the  test  strip. 
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Figure  2.  2-14  Scene  Matching  Area  Correlator 


Since  a  reference  image  must  be  stored  for  each  fix  taken,  it 
is  possible  to  obtain  only  a  few  fixes  with  an  optical  correlation  device. 
Therefore,  the  limiting  accuracy  of  an  optical  correlation/inertial  system 
can  be  considered  to  be  the  position  fix  accuracy  of  the  optical  device. 

2. 2. 9  Doppler/lnertial  Guidance 

.  Unlike  the  other  guidance  systems  discussed  in  this  document, 
the  doppler/inertial  system  does  not  provide  position  updates.  It  was 
assumed  in  the  present  study  that  the  target  location  was  known  without 
error  at  the  beginning  of  the  flight  and  therefore  the  limiting  accuracy 
is  simply  the  position  error  growth  during  the  flight.  The  limiting  veloc  - 
ity  error  for  the  doppler/inertial  system  is  composed  of  a  bias  and  scale 
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factor  error  in  the  doppler  measurement.  (Uncorrelated  measurement 
noise  can  be  averaged  out  if  many  measurements  are  taken. )  The  position 
error  is  the  integral  of  the  velocity  measurement  error  over  the  time  of 
flight  as  expressed  by  Eq.  (2.  2-14). 

Tf 

«P.=  f  (5Vb  +  Vm6SF)dt  (2.2-14) 


where 

6P  =  position  error 

fiV^  =  bias  error  in  doppler  measurement 

V  =  missile  velocity 
m 

8SF  =  doppler  measurement  scale  factor  error 
»  missile  flight  time 


If  the  bias  and  scale  factor  errors  in  the  doppler  measurement 
are  assumed  uncorrelated  the  rms  value  of  the  position  error  becomes 


(2.  2-15) 


where 

6^  =  rms  position  error 
6y^  =  rms  bias  error  in  doppler  measurement 
dgj,  =  rms  doppler  measurement  scale  factor  error 

For  example,  if  the  doppler  bias  error  is  0. 1  ft/sec,  the  doppler  scale 
factor  error  is  10~^,  the  missile  velocity  is  2,000  ft/sec,  and  the  mis¬ 
sile  flies  for  200  sec,  evaluation  of  Eq.  (2.  2-15)  indicates  an  rms  value 
of  the  limiting  position  error  due  to  doppler/inertial  errors  of  about  45  ft. 
As  noted  previously  the  determination  of  the  limiting  accuracy  for  the 
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doppler/inertial  system  was  based  on  the  assumption  of  perfect  knowl¬ 
edge  of  target  location  at  the  beginning  of  the  flight.  Since  the  time  of 
flight  of  a  tactical  missile  is  typically  short,  correlation  between  position 
and  velocity  errors  can  be  ignored.  Therefore,  the  total  delivery  error 
may  be  considered  as  the  root  sum  square  (rss)  of  the  initial  target  loca¬ 
tion  error  and  the  position  error  growth  during  the  flight.  * 

2. 3  CHAPTER  SUMMARY 

This  chapter  has  presented  the  analysis  of  the  limiting  ac¬ 
curacies  for  nine  aided  inertial  guidance  techniques.  The  factors  shov/n 
to  limit  each  technique  are  summarized  in  Table  2. 3-1.  The  effects  of 
these  factors  on  each  technique  are  discussed  individually  in  this  chapter. 


TABLE  2.  3-1 

FACTORS  LIMITING  GUIDANCE  ACCURACY 


Geophysical 

Uncertainties 

Geometric 

Dilution 

Propagation 

Uncertainties 

Electronic 
Clock  . 
Accuracy 

Other 

External  Aid 
Fix  Errors 

Pure  Inertial 

X 

LORAN/Inertlal 

X 

X 

Direct  Ranging 
LORAN/Inerttal 

X 

X 

X 

OMEGA/ Inertial 

X 

X 

DME/Incrtial 

X 

X 

Satellite/Inertial 

X 

X 

X 

Radar  Correlation/ 
beitial 

X 

Optical  Correlation/ 
Inertial 

X 

Doppler/inertial 

X 

A  ' 

This  technique  is  valid  only  if  the  initial  target  location  error  is  not 
dependent  on  the  initial  missile  velocity  error. 
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3.  ACCURACY  SENSITIVITY  STUDIES  FOR 

AUTONOMOUS  MISSILE  GUIDANCE 


3. 1  INTRODUCTION 

This  chapter  discusses  the  accuracy  sensitivity  analysis  of 
three  inertial  guidance  techniques  for  air-to-surface  missiles.  These 
guidance  techniques  are  position  reference  guidance  (inertial  navigation), 
heading  command  guidance,  and  wind  compensated  heading  command 
guidance.  These  analyses  relate  the  sensitivity  of  missile  guidance  er¬ 
rors  to  sensor  errors  and  initial  alignment  errors  for  each  technique. 

Five  representative  trajectories  were  considered  in  the  analyses.  These 
trajectories  were  defined  to  provide  typical  missile  acceleration  pro¬ 
files  so  that  the  effect  of  these  accelerations  on  guidance  errors  could 
be  studied. 

The  implementations  assumed  for  the  three  inertial  guidance 
techniques  are  as  follows: 

Position  Reference  —  The  inertial  navigator  considered  is 
representative  of  a  north,  east,  vertical  inertial  platform  with  two  ac¬ 
celerometers  and  three  gyros.  A  vertical  channel  for  the  inertial  system 
was  not  considered.  It  was  assumed  that  altitude  information  was  de¬ 
rived  from  a  separate  source  such  as  a  barometric  altimeter. 

Heading  Command  Guidance  —  Heading  command  guidance  is 
considered  to  be  performed  by  commanding  predetermined  missile  altitudes 
and  headings  with  respect  to  a  platform  stabilized  in  inertial  space. 
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(In  this  analysis  altitude  information  is  assumed  to  be  derived  from  an 
external  source. )  The  platform  is  assumed  to  be  gimballed  with  three 
gyros  mounted  on  it  to  provide  stabilization.  The  input  axes  of  these 
gyros  are  considered  to  be  mutually  orthogonal. 

Wind  Compensated  Heading  Command  Guidance  —  Wind  com¬ 
pensated  heading  and  altitude  command  guidance  is  similar  to  heading 
command  guidance  discussed  above;  It  differs  only  in  that  the  magnitude 
of  the  wind  velocity  is  measured  prior  to  missile  launch  and  this  is  com¬ 
pensated  for  throughout  the  missile  flight.  The  effect  of  wind  on  missile 
accuracy  is  dependent  on  the  quality  of  the  wind  measurement  and  the 
spatial  variation  in  magnitude  of  the  wind. 

The  following  sections  present  a  description  of  the  trajectories 
considered,  a  discussion  of  inertial  sensor  errors,  and  the  results  of  the 
analyses  of  the  guidance  techniques. 

3. 2  TRAJECTORIES 

Five  representative  trajectories  which  might  be  flown  by  an 
air-to-surface  missile  were  considered  in  the  sensitivity  analyses. 

These  are  discussed  separately  in  the  following  text. 

3. 2. 1  Cruise-GLMe  Trajectory 

The  cruise -glide  trajectory  is  illustrated  in  Fig.  3. 2-1.  It 
consists  of  a  152  sec  straight  and  level  cruise  flight  at  30, 000  ft  altitude, 
followed  by  a  1-g  tip-over  maneuver  to  a  pitch  angle  of  -10  deg.  The 
missile  then  glides  at  this  attitude  until  it  reaches  an  altitude  of  5,000  ft. 
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CONSTANT  MISSILE  STEED  (2000  (</■«') 
CONSTANT  RADIUS  la  TU*-CV£R  (0016iod/i«) 


Not*'-  Eigix*  n«t  10  Koto 


TIME  »0»«  >52  too 


Figure  3. 2 -1  Cruise-Glide  Trajectory 


The  missile  speed  is  constant  (2,000  ft/sec)  and  it  is  heading  north* 
during  the  entire  trajectory.  Examination  of  the  trajectory  shows  that 
the  missile  accelerates  only  during  the  tip-over  maneuver.  Therefore, 
acceleration  dependent  errors  are  only  excited  during  this  maneuver. 


3. 2. 2  Glide -Cruise -Popup  Trajectory 

The  glide-cruise  popup  trajectory  is  shown  in  Fig.  3.  2-2. 

It  consists  of  a  80  sec  glide  at  a  pitch  angle  of  -10  deg  during  which  the 
missile  descends  from  30,000  ft  altitude  to  about  500  ft.  This  is  followed 
by  a  1-g  pull-up  maneuver  to  level  and  a  78  sec  straight  and  level  cruise 
flight.  The  missile  then  performs  a  popup  maneuver  in  which  it  does  a 
4-g  pull-up  to  45  deg  followed  by  a  2-g  tip-over  to  return  the  missile  to 


The  selection  of  the  north  heading  in  this  trajectory  and  others  which 
follow  is  only  for  convenience  and  other  headings  will  not  significantly 
alter  the  magnitude  of  the  results,  only  the  direction  of  the  errors. 


3-3 


CONSTANT  Mi&SlU  iHt 0  12000  It/**) 


Nou  t  figwf*  imm  to  uoU 


*•$36$ 


CONSTANT  RADIUS  ?g 


Figure  3. 2-2  Glide-Cruise-Popup  Trajectory 

an  altitude  of  500  ft.  During  the  tip -over  the  missile  pitch  attitude  changes 
101  deg.  The  missile  speed  is  constant  (2,000  ft/sec)  and  it  is  heading 
north  during  the  entire  trajectory.  This  trajectory  contains  several  high 
acceleration  pvll-up  and  tip-over  maneuvers,  causing  acceleration  depend 
errors  to  be  excited  during  much  of  the  trajectory. 

3. 2. 3  Loft -Cruise -Glide  Trajectory 

The  loft -cruise -glide  trajectory  is  illustrated  in  Fig.  3. 2-3. 

It  consists  of  a  5  sec  constant  acceleration  boost  phase  at  a  missile  pitch 
angle  of  45  deg.  During  this  boost  the  missile  velocity  increases  from 
1,000  ft/sec  to  4,000  ft/sec.  This  is  followed  by  a  loft  phase  during 
which  no  missile  thrust  is  applied  and  the  only  forces  cn  the  missile  are 
due  to  drag,  gravity,  and  lift.  A  drag  model  was  chosen  in  which  the 
dynamic  pressure  and  therefore  drag  force  was  proportional  to  the  square 
of  the  missile  velocity  (Ref.  8).  Using  this  model  an  acceleration  profile 
was  obtained  which  would  deliver  the  missile  to  an  approximate  altitude 
of  25,000  ft  with  a  horizontal  velocity  of  2,000  it/sec.  The  acceleration 
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Figure  3.  2-3  Loft -Cruise -Glide  Trajectory 


profile  obtained  is  shown  in  Fig.  3. 2-4.  *  After  the  loft  phase  the  missile 
maintains  straight  and  level  flight  for  147  sec  follow  by  a  1-g  tip-over 
maneuver  to  a  pitch  angle  of  -10  deg.  The  missile  then  glides  at  this 
angle  to  an  altitude  of  5,000  ft.  After  the  loft  phase  the  missile  velocity 
is  constant  (2,000  ft/sec)  and  during  the  entire  trajectory  the  missile 
is  heading  north. 


The  actual  acceleraions  experienced  by  a  missile  are  a  function  of  many 
parameters  such  as  missile  surface  area,  drag  coefficient,  guidance  law, 
etc.  This  profile  was  chosen  since  it  was  felt  to  be  typical  of  a  realistic 
missile  acceleration  profile.  The  purpose  of  including  missile  accelera¬ 
tions  in  this  analysis  was  to  provide  a  realistic  representation  of  the  ac¬ 
celerations  which  excite  ceitain  sensor  errors.  It  is  felt  that  other  ac¬ 
celeration  profiles  will  excite  these  sensor  errors  in  a  similar  manner 
and  that  they  would  cause  no  significant  deviation  in  the  end  result. 
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Figure  3.  2-4  Acceleration  Profile  During  Loft  Phase  of 
Loft-Cruise-Glide  Trajectory 

> 


3.  2.  4  Loft -Glide-Cruise  Trajectory 


The  loft -glide -cruise  trajectory  is  illustrated  in  Fig.  3.2-5. 

It  consists  of  a  boost  and  loft  phase  similar  to  that  discussed  for  the  loft- 
cruise-glide  trajectory.  This  is  followed  by  a  1-g  tip-over  maneuver  to 
a  pitch  angle  of  -10  deg.  The  missile  then  glides  for  60  seconds  and  per¬ 
forms  a  1-g  pull-up  to  level.  It  then  maintains  straight  and  level  flight 
at  an  altitude  of  500  ft  for  129  sec.  After  the  loft  phase,  the  missile 
velocity  is  constant  (2,000  ft/sec)  and  during  the  entire  trajectory  the 
missile  is  heading  north.  ^ 


3. 2.  5  Turn  After  Launch  With  Terminal  Popup  Trajectory 

A  trajectory  was  studied  in  which  the  missile  performed  a 
turn  after  launch  with  a  popup  maneuver  at  the  end  of  the  trajectory. 
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Figure  3.  2-5  Loft -Glide-Cruise  Trajectory 

This  trajectory  is  illustrated  in  Fig.  3.2-6.  The  missile  is  initially 
heading  north  but  .  after  about  91  sec  of  straight  flight  it  makes  a  1-g 
45  deg  turn  to  the  left.  The  terminal  popup  maneuver  was  similar  to 
the  one  discussed  in  Section  3.  2-2  for  the  glide-cruise-popup  trajectory. 
It  consisted  of  a  4-g  pull-up  to  45  deg  followed  by  a  2-g  tip-over  to  return 
the  missile  to  an  altitude  of  1, 000  ft.  During  the  entire  trajectory  the 
missile  speed  was  constant  (2,000  ft/sec). 

CONSTANT  MISSILE  SPEED  IJOOO  lt/i«)  Note  Flour,  not  to  icol. 

TOTAL  HORIZONTAL  DISTANCE  TRAVELED  *  75nm 

CONSTANT  RADIUS  20 
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3.  3  INERTIAL  SENSOR  ERRORS 


Gyros  and  accelerometers  are  not  perfect  instruments.  Errors 
in  these  instruments  may  contribute  significantly  to  position,  velocity 
and  attitude'  errors,  especially  when  missile  maneuvers  excite  the  ac¬ 
celeration-sensitive  terms.  In  the  sensitivity  analyses  which  were  per¬ 
formed,  inertial  sensor  errors  were  grouped  according  to  type  and 
the  effect  of  each  group  was  then  investigated  separately.  For  example, 
the  influence  of  all  the  mass  unbalance  gyro  errors  was  treated  as  a 
single  effect,  rather  than  investigating  mass  unbalance  errors  in  the 
east  gyro,  etc.  What  follows  is  a  list  of  inertial  sensor  errors  con¬ 
sidered  and  an  explanation  of  the  grouping  used  in  error  sensitivity. cal¬ 
culations.  The  accelerometer  errors  discussed  apply  only  to  position 
reference  guidance  since  heading  command  guidance  requires  only  gyros. 


Accelerometer  Bias  Errors  —  Each  of  the  two  inertial  navi¬ 
gator  accelerometers  is  considered  to  have  a  bias  error  associated  with 
it.  These  errors  contribute  to  velocity  errors  throughout  the  entire 
flight. 

Accelerometer  Scale  Factoif  Errors  —  These  errors  generate 


errors  in  the  sensed  acceleration  proporational  to  the  specific  force 
components  which  lie  along  the  instruments'  input  axes.  Errors  do  not 
usually  appear  during  unaccelerating  flight,  but  this  error  source  may 
become  significant  during  high  acceleration  maneuvers  such  as  turns  or 
a  missile  boost  phase. 


Accelerometer  Nonlineariti! 


errors  in  sensed  acceleration  that  ar^  proportional  to  the  square  and 


ies  —  These  imperfections  create 
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cube  of  the  accelerometer  inputs.  As  with  scale  factor  errors  there 
effects  appear  during  missile  maneuvers. 

Accelerometer  Misalignments  —  This  group  consists  of  two 
effects  due  to  misalignment  of  the  inertial  navigator  accelerometer  input 
axes  —  one  for  each  of  the  lerel  accelerometers.  In  the  truth  model 
two  accelerometer  misalignment  terms  are  omitted  because  the  input 
axes  of  the  level  accelerometers  define  the  navigator  platform  attitude 
achieved  during  alignment.  Errors  in  the  sensed  acceleration  are  gen¬ 
erated  by  misresolution  of  horizontal  accelerations  through  the  misalign¬ 
ment  angle,  causing  this  effect  to  appear  only  during  missile  manuevers 
which  generate  horizontal  accelerations. 

Gyro  Drift  Bias  Errors  —  Each  of  the  three  gyro  input  axes* 
is  considered  to  have  a  bias  drift  rate  error  associated  with  it.  These 
errors  (especially  the  azimuth  gyro  drift  rate)  can  cause  significant 
attitude  errors.  The  attitude  errors  cause  misresolution  of  the  specific 
force  vector,  generating  velocity  and  position  errors. 

Gyro  Scale  Factor  Errors  —  Each  gyro  input  axis  is  considered 
to  have  a  drift  rate  error  proportional  to  the  command  platform  angular 
rate  about  that  axis.  This  generally  results  from  imperfect  calibration 
of  the  gyro  torquers.  There  is  one  scale  factor  error  term  for  each  of 
the  three  gyros  in  the  platform. 

Gyro  Mass  Unbalance  Effects  —  Gyro  drift  rates  proporational 
to  specific  force  components  which  lie  in  the  direction  of  the  spin  and  input 
axes  can  be  explained  in  terms  of  mass  unbalance  effects.  The  salient 

* 

Three  gyros  a^e  assumed  on  both  the  inertial  navigator  and  stable  plat¬ 
form  (used  for  heading  command  guidance). 
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characteristic  of  this  drift  rate  term  is  that  it  is  directly  proportional 
to  specific  force.  Since  preflight  alignment  and  calibration  will  permit 
compensation  of  the  effects  of  a  1-g  environment,  mass  unbalance  drift 
rate  terms  will  contribute  to  guidance  errors  only  during  missile  maneu¬ 
vers. 

Gyro  Anisoelasticity  Effects  —  Gyro  drift  rates  proportional 
to  the  second  power  of  acceleration  result  from  these  effects.  Major 
anisoelastic  effects  result  from  simultaneous  specific  force  components 
along  the  spin  and  input  axes.  They  are  excited  only  during  missile 
maneuvers,  such  as  turns  and  pullups. 

3.  4  POSITION  REFERENCE  (INERTIAL)  GUIDANCE 

The  inertial  navigator  considered  in  this  analysis  is  repre¬ 
sentative  of  a  north,  east,  down  inertial  platform  with  two  accelerometers 
and  three  gyros.  A  vertical  channel  for  the  inertial  system  is  not  con¬ 
sidered;  altitude  information  is  assumed  to  be  derived  from  a  separate 
source.  This  section  presents  the  inertial  sensor  error  sensitivity 
sensitivity  tables  for  the  position  reference  system  for  the  five  trajec¬ 
tories  discussed  previously.  The  use  of  the  tables  provides  only  the 
guidance  error  sensitivity  to  inertial  sensor  errors  and  platform  align¬ 
ment  errors.  Other  sources  of  error  will  also  be  present  for  the  posi¬ 
tion  reference  technique,  such  as  the  effect  of  vertical  deflections  and 
gravity  anomalies  on  the  inertial  navigator.  These  errors  have  not  been 
explicitly  discussed  in  this  section  but  their  influence  may  be  calculated 
by  use  of  Fig.  2. 1-5  in  Section  2. 1-3. 
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3. 4. 1  Inertial  Sensor  Error  Sensitivity  Tables 

The  sensitivity  tables  for  the  guidance  errors  ai  the  end  of  the 
five  trajectories  are  presented  in  this  section.  Table  3.  4-1  presents  a 
summary  of  the  trajectory  and  corresponding  table  in  which  the  results 
may  be  found.  The  sensitivity  tables  give  the  north  and  east  position 
and  velocity  errors  along  with  the  CEP  of  the  position  errors.  These 
results  are  given  for  unit  sensor  error  inputs  and  in  this  form  are  not 
representative  of  a  particular  navigation  unit.  To  use  the  results  the 
values  can  be  scaled  by  the  proper  sensor  error  value,  due  to  the  lin¬ 
earity  of  the  error  equations.  An  example  of  this  scaling  is  presented 
in  the  following  section.-.. 


TABLE  3. 4-1 


INDEX  FOR  SENSITIVITY  TABLES 


Trajectory 

Table  Number 
of  Sensitivity  Table 

Cruise- Glide 

3. 4-2 

Glide-Cruise- Popup 

3.4-3 

Loft- Cruise-Glide 

3.  4-4 

Loft- Glide- Cruise 

3.  4-5 

4  • 

Turn  After  ^aunch 

3.4-6 

with  Terminal  Popup 

Minus  sign  indicates  values  to  be  squared  and  subtracted  when  rss  of  error  sources  is  performed  to  determine  the  total  error 


Minas  sign  indicates  values  to  be  squared  and  subtracted  “'hen  rss  of  error  sources  is  performed  to  determine  the  total  error 


.  ■■  * 


Minus  cign  indicates  values  to  be  squared  and  subtracted  when  rss  of  error  sources  is  performed  to  determine  the  total  error 
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3. 4. 2  Sample  Calculation 

The  results  of  the  sensitivity  analysis  can  best  be  used  by  refer¬ 
ring  to  a  particular  navigation  system.  In  this  section,  a  typical  missile 
navigation  system  is  described  and  the  guidance  accuracies  achieved  with 
the  previously  discussed  trajectories  are  computed.  These  calculations 
demonstrate  the  use  of  the  sensitivity  tables  and  some  comparisons  which 
can  be  made  using  the  results. 

A  particular  alignment  scheme  and  a  navigation  system  which 
would  by  typical  of  a  low  cost  navigator  for  use  in  a  missile  application 
were  assumed.  The  error  coefficients  assumed  are  given  in  Table  3. 4-7. 

TABLE  3.  4-7 

ALIGNMENT  AND  NAVIGATOR  ERRORS 


Error  Source 

Error  Magnitude 

Initial  position  errors 

200  ft 

Initial  velocity  errors 

2  ft/sec 

Initial  tilt  errors 

0.1rar 

Initial  azimuth  error 

1:0  rar 

Accelerometer  bias  errors 

5  x  10‘4  g 

Accelerometer  bias  error  correlations 
with  initial  tilt  error 

0.  9  correlation  coef. 

Accelerometer  scale  factor  errors 

500  ppm  . 

Accelerometer  misalignment  errors 

1  mr 

Accelerometer  nonlinearities  (g  squared) 

8  x  10'5  g/g2 

* 

•Accelerometer  nonlinearities  (g  cubed) 

5  x  10‘6  g/g3 

Gyro  drift  bias  errors 

1  deg/hr 

Gyro  drift  bias  error  correlations 
with  initial  azimuth  error 

0. 9  correlation  coef. 

Gyro  scale  factor  errors 

1000  ppm 

Gyro  mass  unbalance  errors 

3. 0  deg/hr/g 

Gyro  anisoelasticity 

0. 3  deg/hr/g2 

■* 

Values  not  presently  available  for  LCP  HI  platform,  typical  values 
vere  assumed 
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Except  in  the  cases  noted  the  sensor  performance  data  for  the  Raytheon 
LCP-m  low  cost  inertial  platform  was  used.  The  alignment  scheme  to 
which  the  numbers  correspond  might  be  a  transfer  alignment  from  a  high 
quality  aircraft  navigator  with  a  position  update  prior  to  laupch. 

The  guidance  error  CEP's  for  the  five  trajectories,  based  on 
the  assumed  error  magnitudes  are  presented  in  Table  3.  4-8.  These  were 
computed  by  scaling  the  sensitivity  results  by  the  error  magnitudes  as 
demonstrated  in  the  following  text.  Using  Table  3. 4-2  which  contains  the 
sensitivity  results  for  the  cruise-glide  trajectory,  the  following  calcu¬ 
lations  can  be  performed.  From  the  table  it  can  be  seen  that  a  1  ft/sec 
initial  velocity  error  causes  a  226  ft  position  CEP.  The  initial  velocity 
error  for  the  alignment  scheme  assumed  was  2  ft/sec  (Table  3.  4-7). 

The  resultant  delivery  error  CEP  is  therefore 

[(266  ft)  /  (1  ft/sec)]  (2  ft/sec)  =  532  ft 

All  other  sensitivity  calculations  are  performed  in  a  similar  manner  with 
one  exception.  When  scaling  the  terms  representing  correlations  between 
accelerometer  bias  error  and  initial  tilt  error  or  gyro  drift  bias  error 
and  initial  azimuth  error  the  following  method  must  be  used:  Assume 
b  is  the  resultant  position  or  velocity  error  caused  by  a  correlation  be¬ 
tween  error  sources  Xj  and  Xg.  If  the  rms  magnitude  of  x^  or  x9  or  the 
correlation  coefficient  for  tire  guidance  system  being  analyzed  differs 
from  the  unit  input  value  then  the  result,  b,  must  be  scaled  as  follows.  * 

-  *  ~ 

This  scaling  results  since  correlation  terms  are  off-diagonal  elements 
in  the  covariance  matrix. 
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TABLE  3.  4-8 

CEP  OF  GITOANCE  POSITION  ERRORS  FOR  5  TRAJECTORIES 


Error  Source* 

CEP 

Cruise 

tilde 

(10 

CEP 

Clide 

Cruise 

Popup 

(ft) 

CEP 

Loft 

Cruise 

Glide 

(ft) 

CEP 

Loft 

Glide 

Cruise 

(ft) 

CLP 

Turn  After 
Launch  with 
Terminal 
Popup 

(ft) 

Initial  position  errors 

236 

236 

232 

232 

235 

Initial  velocity  errors 

S3  2 

550 

030 

532 

548 

Initial  tilt  errors 

96 

110 

78 

78 

102 

Initial  azimuth  error 

3 

6 

160 

169 

105 

Accelerometer  bias  error* 

464 

530 

402 

493 

.525, 

Accelerometer  bias  error  correlation* 

-206 

-220 

-186 

-185 

-220 

Accelerometer  scale  factor  errors 

,1 

M 

m 

83 

S3 

68 

Accelerometer  misalignment  error* 

1 

4 

166 

166 

136 

Accelerometer  nonllnea-ltles  (g  squared) 

— 

•  1 

297 

297 

9 

Accelerometer  nonltnearittes  (7  cubcdi 

— 

— 

229 

229 

— 

Cyro  drift. bias  errors 

344 

406 

315 

351 

331 

Gyro  drift  bias  error  correlation* 

•  13 

-11 

-12 

-12 

-122 

Gyro  scale  factor  errors 

6 

7 

5 

0 

■  6 

Gyro  mass  unbalance  errort 

1 

6 

297 

292 

178 

Gyro  anisoclastlrity 

— 

— 

34 

35 

5 

rrs  Total 

817 

873 

969 

983 

890 

Minus  sign  Indicates  values  that  were  squared  and  subtracted  when  rss  total  was  computed 


—  Indicates  error  sources  which  had  negligible  contribution 


! 

i 


1 -  i  /  rmi  value  of  x^ 

/  rms  value  of 

actual  correlation 
coefficient 

“  ^aiue*0  M / luUl  inl5ut  value  °*  xl 

J  unit  input  value  of  x ^ 

input  correlation 
coefficient 

(3.  4-1) 


To  clarify  this  scaling  calculation  consider  the  following  example:  From 
Table  3. 4-2  for  the  cruise-glide  trajectory,  the  correction  to  position 
CEP  caused  by  a  0. 9  correlation  between  a  1-g  accelerometer  bias 


-  3-19 

I 


THE  ANALYTIC  SCIENCES  CORPORATION 


4 

error  and  a  1  mr  initial  tilt  error  is  -2. 92  x  10  ft.  *  For  the  example 
considered  in  Table  3.  4-7,  the  correlation  coefficient  ir.  0.9,  the  ,>c~ 
celeror.ieter  bias  error  is  5  x  10"^  g,  and  the  initial  tilt  error  ir  0. 1  mr. 
Using  the  scaling  in  expression  (3.  4-1),  the  resultant  correction  lo  posi¬ 
tion  CEP  for  the  example  being  considered  is 


(-2.  92  x  10*  ft) 


'0.  lmr 
I.  Omr 


5  x  10"4  g  (0.  9) 


206  ft 


What  follows  is  a  partial  list  of  conclusions  which  can  be  ob¬ 
tained  from  these  results,  to  demonstrate  their  usefulness.  It  should  be 
noted  that  the  conclusions  reached  do  not  apply  in  general  but  arc  depen - 
ent  on  the  input  val'  es  assumed  in  this  example. 


•  Initial  velocity  errors  arc  the  largest  single 
error  contributor  for  all  trajectories.  •  Use  of 
improved  initial  velocity  information  (this  could 
be  accomplished  by  doppler  aiding  the  master 
aircraft  navigator  before  transfer  alignment) 
could  significantly  reduce  this  error  source. 

For  example  if  the  initial  velocity  error  were 
reduced  to  0. 2  ft/sec  then  the  contribution  of 
velocity  error  would  b?  reduced  by  a  factor 

of  10  in  each  ease. 

•  For  all  trajectories  at  least  85%  of  the  total 
rss  error  was  due  to  initial  alignment  errors, 
accelerometer  bias  error,  and  gyro  diift  errors. 
Improved  accuracy  could  be  obtained  by  improving 
the  alignment  scheme  and/or  using  inertial  sensors 
with  smaller  bias  errors. 


The  minus  sign  indicates  that  the  value  is  to  be  squared  and  subtracted 
when  the  total  rss  error  is  computed.  The  minus  sign  shows  that  the  cor¬ 
relation  between  two  error  sources  helps  reduce  their  total  effect. 
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•  In  the  two  trajectories  containing  a  loft  phase, 
the  acceleration  dependent  errors  —  acceler¬ 
ometer  nonlinearities,  gyro  mass  unbalance, 
etc.  —  generate  large  guidance  errors.  Due 
to  the  higher  contribution  of  other  error  sources 
discussed  previously  these  acceleration  depen¬ 
dent  errors  do  not  significantly  degrade  overall 
performance  when  compared  to  trajectories 
which  contain  few  high  accelerations.  There¬ 
fore,  use  of  improved  inertial  sensors  which 
have  smaller  error  coefficients  for  acceler¬ 
ation  dependent  terms  would  provide  no  signi¬ 
ficant  guidance  accuracy  improvement. 

The  sensitivity  tables  presented  in  this  section  provide  a 
simple  tool  for  performing  accuracy  tradeoff  studies.  They  allow  the 
evaluation  and  comparison  of  mar'  systems  in  a  short  period  of  time 
and  give  insights  into  areas  in  which  component  or  procedural  changes 
would  be  most  significant. 


3.  5  HEADING  COMMAND  GUIDANCE 

This  section  discusses  the  sensitivity  analysis  of  the  heading 
command  and  wind  compensated  heading  command  guidance  techniques. 
These  techniques  differ  only  in  the  manner  in  which  wind  effects  the 
guidance  errors;  all  other  error  contributors  produce  the  same  guid¬ 
ance  errors  for  the  two  techniques.  Both  guidance  techniques  are  con¬ 
sidered  to  be  implemented  by  commanding  predetermined  missile  altitudes 
and  headings  with  respect  to  a  platform  stabilized  in  inertial  space. 

The  platform  is  assumed  to  be  gimballed  with  three  single -degree-of- 
freedom  gyros  mounted  on  it  to  provide  stabilization  as  illustrated  in 
Fig.  3.5-1. 
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H-93SS 

SA-SPIN  AXIS 
JA- INPUT  AXIS 


Figure  3.  5-1  Gyro  Configuration  Assumed  for  Heading 
Command  Guidance 


The  input  axies  of  these  gyros  are  considered  to  be  mutually  orthogonal.  * 
Altitude  information  is  assumed  to  be  derived  from  an  external  source. 

The  inertial  sensor  error  sensitivity  tables  for  the  heading 
command  guidance  techniques  are  presented  for  the  five  trajectories 
discussed  previously.  The  use  of  the  tables  provides  only  the  guidance 
error  sensitivity  to  inertial  sensor  errors  and  platform  alignment  errors. 
Another  source  of  error  important  for  the  heading  command  guidance 
techniques,  the  effect  of  wind,  is  presented  in  a  separate  part  of  this 
section.  Finally  a  sample  calculation  is  presented,  demonstrating 
the  use  of  the  analyses  to  evaluate  a  particular  system. 


* 

A  second  gyro  configuration  may  be  used.  For  this  configuration  the 
spin  axis  of  the  east  gyro  is  vertical.  It  was  found  that  this  change  in 
configuration  does  not  significantly  change  the  guidance  errors  for  the 
five  trajectories  considered. 
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3.  5. 1  Inertial  Sensor  Error  Sensitivity  Tables 

The  sensitivity  tables  fou  the  heading  command  guidance  errors 
at  the  end  of  the  five  trajectories  are  presented  in  this  section.  Table  3*  5-1 
presents  a  summary  of  the  trajectories  and  corresponding  tables  in  which 
the  results  may  be  found.  The  sensitivity  tables  give  the  along  and  across 
track  position  errors  along  with  the  CEP  of  the  pos  tion  errors. 

These  results  are  given  for  unit  sensor  error  inputs  and  in 
this  form  are  not  representative  of  a  particular  navigation  unit.  To  use 
the  results  the  values  can  be  scaled  by  the  proper  sensor  error  value 
(as  illustrated  in  Section  3.  4. 2),  due  to  the  linearity  of  the  error  equa¬ 
tions. 

For  convenience,  one  error  contributor  which  is  not  attributed 
to  inertial  sensors  has  been  included  in  the  sensitivity  tables.  This  error 
source  is  the  error  in  the  knowledge  of  the  missile  velocity  and  would 
results  from  errors  in  the  missile  speed  indicator  or  thrust  control. 

For  a  heading  command  guidance  system  which  does  not  rely  on  external 
position  or  velocity  information,  the  along  track  missile  position  must 
be  determined  from  a  knowledge  of  the  missile  velocity  and  time  of  flight. 

If  there  is  an  error  in  the  knowledge  of  the  missile  velocity,  an  error 
in  the  along  track  position  will  be  generated.  When  external  position  or 
velocity  information  is  available,  this  error  will  not  be  significant.  In 
the  present  analyses  the  missile  velocity  error  is  considered  to  be  directly 
proportional  to  the  missile  velocity. 

3.  5. 2  Wind  Compensr  don 

One  of  the  major  sources  of  guidance  error  when  attitude  com¬ 
mand  guidance  is  used  in  a  cruise  missile  is  the  effect  of  wind.  It  has 
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TABLE  3.  5-1 

INDEX  FOR  SENSITIVITY  TABLES 


Trajectory 

Table  Number 

Cruise -Glide 

3.  5-2 

Glide-Cruise-P  opup 

3.5-3 

Loft -Cruise -Glide 

3.  5-4 

Loft -Glide-Cruise 

’  3.5-5 

Turn  After  Launch 

3.  5-6 

with  Terminal  Popup 

been  suggested  that  the  wind  be  somehow  measured  and  the  commanded 
heading  or  velocity  be  compensated  to  remove  the  effects  of  wind.  If  a 
perfect  measurement  of  the  wind  is  achieved  and  if  the  wind  is  constant 
after  the  measurement,  the  influence  of  wind  on  guidance  errors  can  be 
eliminated.  In  truth  the  wind  is  seldom  constant  and  a  perfect  measure¬ 
ment  is  not  possible.  This  section  analyzes  the  effect  of  the  wind  on  mis¬ 
sile  guidance  errors  for  wind  compensated  and  uncompensated  heading 
command  guidance  systems  when  the  wind  is  not  constant  over  the  range 
of  the  missile.  (Temporal  variations  are  ignored  in  favor  of  spatial 
changes  in  wind  velocity. )  The  wind  is  further  assumed  to  be  a  homo¬ 
geneous,  isotropic  process  allowing  the  effect  of  the  wind  on  crosstrack 
and  along  track  guidance  errors  to  be  treated  separately. 
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If  the  wind  velocity*,  v  ,  is  represented  as  a  first-order  markov 

process  in  spatial  coordinates,  it  and  the  resultant  guidance  error  due  to 

wind,  e  ,  are  related  by  the  following  vector -matrix  differential  equation 
w 

where  T  is  the  characteristic  time  of  the  random  process. 


(3.  5-1) 


For  a  missile  with  constant  velocity**,  V,.  T  =  V/D  where  D  is  the  charac¬ 
teristic  distance  of  the  wind  markov  process.  The  variable  u  is  white 
noise  forcing  the  first-order  linear  system  whose  output  is  wind: 


E[u(t)u(r)]  =  6  (t-T) 

2 

where  6  ( )  is  the  Dirac  delta  function  and  cr  is  the  mean  square  wind 

w 

velocity. 

The  covariance  equation  for  the  vector 


* 

The  derivation  which  follows  is  applicable  to  either  an  across  or  along 
track  wind  since  they  may  be  treated  separately. 

*< 

The  effect  of  wind  on  guidance  errors  is  derived  for  a  simplified  case 
of  constant  missile  velocity.  The  results  are  still  generally  applicable 
for  the  five  trajectories  considered  in  this  report  since  the  missile  velocity 
is  constant  throughout  most  of  these  trajectories. 
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in  Eq.  (3. 5-1)  is 

P  =  FP  +  PFT  +  Q  (5.  5-2) 

where  F  is  defined  in  Eq.  (3.  b-1)  and  the  covariance  matrix  P  is  com¬ 
posed  of  elements 

pn  = 

p12  "  E(Vw> 

P«9  =  ^  ^ 

^22  \  w/  =  ew 

0 
0 


and 


Q 


2ct 


w 


From  Eq.  (3.  5-2)  the  differential  equations  for  the  elements  of  p  are 


P11 


2p..  2o  9 

11  w  .  „  /n\  _  2 

— rjT“  +  '  rp  I  Pj  J  (0)  ~  a 


Pi  o 

p12  =  P11  ~  5  P12^0)  =  ° 


(3.  5-3) 


p22  =  2p12 


P22(0)  =  0 
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This  figure  illustrates  the  important  observations  that  follow 
from  Eq,  (3.  5-5). 

It  can  be  seen  that  the  ratio  a/a  is  an  indication  of  bow  ac  - 

w 

curately  the  wind  is  measured  in  a  particular  situation;  when  a/a^  is 

large  a  poor  measurement  of  the  wind  is  av>j  ible.  In  particular 

■  *./; 

a/a  =  1  corresponds  to  no  measurement  O:  ‘we  wind  (uncompensated' 
w  _ 

heading  command  guidance).  Figure  3.  5-2  indicates  that  when  the  mis¬ 
sile  range  is  several  times  the  characteristic  distance  of  the  wind  (b>2) 
there  is  little  value  in  having  a  good  measurement  of  wind.  For  exampl? 
when  b  =  2,  a  is  only  20%  greater  when  no  wind  measurement  is  made 
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y  zy  0.8 


b  v 


Q  .  2  4  6  8  10 

RATIO  Of  RANGE  TO  WIND  CHARACTERISTIC  DISTANCE,  b 


a}  -  MEAN  SQUARE  GUIDANCE  ERROR  DUE  TO  WIND 
ew 

cr*  =  MEAN  SQUARE  WIND 
R  =  RANGE  FROM  TIME  OF  WIND  MEASUREMENT. 


V  =  VELOCITY  (ASSUMED  CONSTANT ) 

a2  *  MEAN  SQUARE  ERROR  IN  MEASUREMENT  OF 

WIND  ( a/aw  =  1  CORRESPONDS  TO  NO  MEASUREMENT 
OF  WIND) 


Figure  3.  5-2  Effect  of  Wind  on  Missile  Guidance  Error  for 
Heading  Command  Guidance 
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than  when  only  10%  error  is  available.  Physically  one  would  expect  the 
value  of  an  initial  measurement  of  wind  to  diminish  as  the  missile  flies 
through  more  and  more  characteristic  distances  of  the  first-order  markov 
process  which  is  used  here  to  represent  wind.  Figure  3.  5-1  reinforces 
this  intuitive  explanation. 

To  illustrate  the  use  of  Fig.  3.  5-2  to  calculate  the  effect  of 
wind  on  the  missile  guidance  errors  the  following  example  is  presented: 
Assume  a  missile  is  flying  with  a  10  ft/stc  wind  acting  on  it  in  the  along 
and  across  track  directions.  If  the  correlation  distance  of  tie  wind  is 
150  nm  and  the  missile  range  is  75  nra  (ihe  missile  range  for  all  fixe 
trajectories  discussed  in  Section  3.  2  is  75  nm),  the  ratio  b  is  1/2. 

For  the  five  trajectories  studied  the  velocity  is  2000  ft/sec  throughout 
most  of  the  trajectory.  Using  Fig.  3.5-2,  the  along  track  or  across 
track  guidance  errors  can  be  found  based  on  the  above  parameters. 

These  are  summarized  in  Table  3.  5-7  along  with  the  position  CEP. 

For  the  particular  example  chos  en  the  wind  correlation  distance  is  suf¬ 
ficiently  large  with  respect  to  the  missile  range  that  significant  guidance 
accuracy  improvement  can  be  achieved  by  wind  compensation. 


3.  5.  3  Sample  Calculation 

As  with  position  reference  guidance,  the  results  of  the  sensi¬ 
tivity  analyses  can  best  be  used  by  referring  to  a  particular  navigation 
system.  In  this  section,  a  typical  missile  navigation  system  is  described 
by  its  assumed  error  coefficients.  The  guidance  accuracies  achieved  with 
the  previously  discussed  trajectories  are  computed.  These  calculations 
demonstrate  the  use  of  the  sensitivity  tables  and  some  comparisons  which 
can  be  made  using  the  results. 
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TABLE  3.  5-7 

EXAMPLE  OF  TII^  EFFECT  OF  WIND 
ON  GUIDANCE  ERRORS 


Guidance  Technique 

Along  and 
Across  Track 
Guidance  Error 
. (ft) . 

- - 

Position 
'  CEP 

(ft) 

Uncompensated  Heading  Command 
,  Guidance 

2100 

2490 

Wind  Compensated  Heading  Command 
Guidance  with  50%  Error  in 

Initial  Wind  Measurement 

1400 

1650 

Wind  Compensated  Heading  Command 
Guidance  with  10%  Error  in 

Initial  Wind  Measurement 

1100  . 

1300 

A  particular  alignment  scheme  and  a  stable  platform  which  would 
be  low  cost  for  use  in  a  missile  application  were  assumed.  The  error  co¬ 
efficients  assumed  were  given  in  Table  3.  5-8.  The  alignment  scheme  to 
which  the  numbers  correspond  might  be  a  transfer  alignment  from  a  high 
quality  aircraft  navigator  with  a  position  update  prior  to  launch.  a 

The  guidance  error  CEP’s  for, the  five  trajectories,  based  on 
the  assumed  error  magnitudes  are  presented  in  Table  3.  5-9.  These  were 
computed  by  scaling  the  sensitivity  results  by  the  error  magnitudes,  as 
demonstrated  in  Section  3. 4. 2.  Table  3.  5-9  includes  only  the  effects  of 
alignment  and  sensor  errors  and  does  not  include  the  effect  of  wind.  This 
.  effect  can  be  computed  separately  based  on  the  particular  wind  parameters, 
as  discussed  in  the  previous  section  and  rss'd  to  give  the  total  guidance 


error. 
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TABLE  3.  5-8 

ALIGNMENT  AND  PLATFORM  ERRORS 


Error  Source 

Initial  position  errors 

Initial  tilt  errors 

Initial  azimuth  error 

Gyro  drift  bias  errors 

Gyro  drift  bias  error  correlations 
with  initial  azimuth  error 

Gyro  scale  factor  errors 

Gyro  mass  unbalance  errors 

Gyro  anisoelasticity 

Velocity  scale  factor  error 


Error  Magnitude 
200  ft 
0.1mr  • 

1. 0  mr 
1  deg/hr 

0. 9  correlation  coef. 

1000  ppm 
3. 0  deg/hr /g 
0. 3  deg/hr /g2 
1000  ppm 


What  follows  is  a  partial  list  of  conclusions  which  can  be  ob¬ 
tained  from  these  results  to  demonstrate  their  usefulness.  It  should  be 
noted  that  the  conclusions  reached  do  not  apply  in  general  but  are  depen¬ 
dent  on  the  input  values  assumed  in  this  example. 

•  0=yro  errors  contribute  less  than  10%  to  the  total 

guidance  error  for  all  trajectories.  Therefore, 
using  more  accurate  gyros  would  provide  no 
significant  improvement. 

•  Comparing  the  results  in  Table  3. 5-9  with  the 
possible  effect  of  wind  (discussed  in  Section  3.  5. 2), 
in  general  it  would  be  expected  that  wind  would 
cause  much  larger  guidance  errors  than  those 

.  contributed  by  platform  sensor  and  alignment 
errors. 
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The  sensitivity  tables  in  this  section  provide  a  simple  tool  for 
performing  accuracy  tradeoff  studies.  They  allow  the  evaluation  and 
comparison  of  heading  command  guidance  systems  in  a  short  period  of 
time.  From  them  insight  can  be  gained  into  areas  In  which  component 
or  precedural  changes  can  most  significantly  reduce  guidance  errors. 
The  reader  is  cautioned  that  when  comparing  heading  command  guidance 
techniques  with  other  guidance  techniques,  the  generally  significant 
contribution  of  wind  must  be  considered  for  the  particular  missile  ap¬ 
plication. 


3. 6  CHAPTER  SUMMARY 

This  chapter  presents  the  accuracy  sensitivity  analysis  for 
three  inertial  guidance  techniques;  position  reference,  heading  command 
guidance,  and  wind  compensated  heading  command  guidance.  Sensitivity 
tables  are  presented  for  five  trajectories  relating  inertial  sensor  errors 
and  initial  alignment  errors  to  missile  guidance  errors.  Sample  calcu¬ 
lations  are  also  presented  to  aid  the  reader  in  using  the  tables.  ’  For  the 
heading  command  guidance  techniques,  a  discussion  of  wind  effects  is 
also  presented  since  this  is  a  potentially  significant  contributor  to  guidance 
errors.  No  detailed  conclusions  as  to  relative  accuracy  of  the  techniques 
can  be  made  since  the  errors  are  dependent  on  the  missile  and  sensor 
parameters.  However,  the  results  in  the  chapter  may  be  used  to  achieve 
detailed  comparison  of  particular  systems. 
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4.  SUMMARY 


4. 1  SUMMARY  AND  CONCLUSIONS 

This  report  presents  the  results  of  two  studies  dealing  with 
the  air -to -surface  missile  guidance  problem.  The  first  portion  covers 
the  study  of  accuracy  limits  for  aided  inertial  guidance  techniques.  This 
involves  the  computation  of  the  lower  bounds  on  the  attainable  accuracy 
of  specific  guidance  techniques.  The  second  part  of  this  document  dis¬ 
cusses  the  sensitivity  studies  which  were  performed  relating  missile 
guidance  errors  to  sensor  and  alignment  errors  for  three  guidance  techni¬ 
ques. 

The  limiting  accuracy  for  nine  proposed  aided  inertial  guid¬ 
ance  techniques  is  determined.  These  systems  cover  a  wide  range  of 
possible  system  designs.  The  systems  studied  are: 

•  Pure  Inertial  Guidance 

•  .  LORAN/lnertial  Guidance 

•  Direct  Ranging  LORAN/Inerlial  Guidance 

•  OMEGA/lndrtial  Guidance  ? 

•  DME/lnertial  Guidance 

•  Sateilite/lnertial  Guidance* 

•  Radar  Correlation/lne.vtial  Guidance* 

Details  are  presented  in  Volume  If  of  this  report,  Ref.  7. 

.  *  ■  _  ‘ 
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•  Optical  Correlation/lnertial  Guidance 

*  ‘  Doppler /inertial  Guidance 

In  all  cases  these  systems  are  assumed  to  employ  an  optimum  infor¬ 
mation  handling  scheme.  This  assumption  provides  a  consistent  base¬ 
line  against  which,  alternative  suboptimal  schemes  can  be  compared  as 
well  as  giving  the  desired  lower  limit  values. 

The  accuracy  of  many  of  the  systems  studied  was  dependent 
on  the  missile  station  geometry.  For  this  reason,  the  results  of  the 
analyses  are  often  presented  in  the  form  of  contour  plots.  These  plots 
provide  contours  of  equal  guidance  error  about  the  station  configuration. 

By  this  method  a  large  amount  of  data  was  conveniently  presented  for  ease 
of  use  by  the  reader. 

The  second  portion  of  this  report  presents  the  sensitivity  studies 
performed  for  three  guidance  techniques;  pure  inertial,  heading  command, 
and  wind  compensated  heading  command  guidance.  The  results  are  pre¬ 
sented  in  the  form  of  sensitivity  tables  which  relate  the  missile  guidance 
errors  to  inertial  sensor  and  alignment  errors  for  five  representative  tra¬ 
jectories.  These  tables  allow  the  rapid  evaluation  and  comparison  of  in¬ 
ertial  sensors  and  alignment  techniques  for  use  in  proposed  guidance  sys¬ 
tems.  To  keep  the  guidance  accuracy  assessment  in  proper  perspective, 
a  discussion  of  another  important  error  source  —  toe  effect  of  wind  -- 
which  is  not  related  to  sensor  or  alignment  errors  has  been  included. 

It  is  difficult  to  draw  conclusions  directly  from  the  sensitivity  tables 
as  presented.--  they  are  best  used  by  referring  to  a  particular  navigation 
system.  For  this  reason  a  sample  calculation  has  been  included  to  de¬ 
monstrate  the  use  of  the  tables  and  some  possible  conclusions  which 
can  be  inferred  from  their  use. 
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4. 2  TECHNOLOGICAL  FORECAST 

This  report  describes  tw.o  valuable  resources  developed  for  the 
analysis  of  missile  guidance  techniques.  The  accuracy  limit  studies  for 
aided  inertial  guidance  techniques  provide  a  lower  limit  for  the  spectrum 
of  guidance  techniques  studied  and  a  consistent  baseline  against  which  al¬ 
ternative  suboptimal  schemes  can  be  compared.  These  analyses  will  aid 
in  the  evaluation  of  proposed  systems  to  establish  the  relative  accuracy 
improvement  they  may  provide.  The  guidance  sensitivity  studies  are 
valuable  in  performing  accuracy  tradeoff  studies  of  inertial  sensors  and 
alignment  techniques  to  establish  possible  guidance  improvements  offered 
by  alternative  instruments.  They  also  provide  an  indication  of  the  areas 
in  which  component  changes  will  cause  the  most  significant  guidance  ac¬ 
curacy  improvements.  These  studies  provide  powerful  reference  material 
for  use  in  evaluating  proposed  guidance  systems,  both  for  the  present  and 
future.  For  this  reason,  summaries  of  these  studies  will  be  incorporated 
into  a  future  version  of  the  tactical  missile  guidance  handbook,  Ref.  9. 
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